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NOMENCLATURE 


Latin  Letters 

Symbol  Definition 

a  The  spherical  area  of  entrance  port  (flux  averaging 
sphere) 

A  Absorptivity 

Aq  .  Absorptivity  with  diffuse  scattering  electrons, 

A^  Absorptivity  with  specular  scattering  electrons. 

A  Bulk  absorptivity. 

B 

b  Effective  prism  base.  Also  spherical  area  of 
detector  port  (flux  averaging  sphere). 

b  Quantum  correction  factor  due  to  electron -phonon 

interaction. 

b°  dc  value  of  quantum  correction  factor  due  to  elec- 
tron-phonon  interaction. 

Bn  Total  flux  reaching  the  detector  port  of  the  flux 
D  averaging  sphere. 

c  Speed  of  light.  Also  spherical  area  of  sample  port 

(flux  averaging  sphere) . 

D  Diameter  of  paraboloid  mirror  in  the  monochromator. 

D,  Recorder  deflection  corresponding  to  SPECIMEN  IN 

position. 

D2  Recorder  deflection  corresponding  to  SPECIMEN  OUT 
position. 

e  Charge  of  an  electron. 

E  Electric  field. 

E^  The  power  emitted  by  the  source  per  unit  area  per 
unit  solid  angle  per  unit  wavelength. 


Definition 


The  power  emitted  by  the  source  per  unit  area  per 
unit  solid  angle  in  the  A  X  range. 

Electron  distribution  function.  Also,  focal  length. 

Planck's  constant/2  it . 

Magnetic  field 

Intensity  after  one  reflection. 

Intensity  after  three  reflections. 

Current  density. 

Absorption  coefficient. 

Boltzmann  constant 

Height  of  the  slit. 

Mass  of  an  electron. 

Index  of  refraction. 

Complex  index  of  refraction. 

The  number  of  electrons  per  unit  volume. 

Fraction  of  electrons  which  are  reflected  specu¬ 
larly. 

Total  radiant  flux  incident  into  the  incident  port 
of  the  flux  averaging  sphere. 

Reflectivity.  Also,  radius  of  the  sphere. 

Reflectivity  in  each  of  the  four  orientations. 

Uncertainty  of  reflectivity. 

Reflectivity  of  standard  mirror. 

Surface  resistance. 

Reflectance  of  the  sample. 

Reflectance  of  sphere  wall. 

Poynting  vector. 


X 


Symbol 

s 

T 

A  T 

W 


V 


v 


f 


X 


X 


z 

z 


ee 


ep 

,° 

ep 


6 

6, 

0 

0. 


D 


Definition 

% 

Mechanical  slit  width. 

Temperature . 

Uncertainty  of  transmission. 

Power  expenditure. 

Electron  velocity. 

Fermi  velocity. 

Surface  reactance. 

Slit  width. 

Penetration  distance  of  the  wave. 

Impedance  of  surface. 

Greek  Letters 

Prism  apex  angle.  Also  polarizability. 

Damping  coefficient. 

Damping  coefficient  due  to  electron-electron  inter 
action . 

Damping  coefficient  due  to  electron-phonon  inter¬ 
action. 

dc  value  of  damping  coefficient  due  to  electron- 
phonon  interaction. 

Real  dielectric  constant. 

Complex  dielectric  constant. 

Depth  of  penetration. 

Skin  depth. 

Angle c 

Debye  temperature. 

Wavelength  {in  microns). 


Definition 


Spectral  slit  width. 

Spectral  slit  width  due  to  dispersion  of  the  prism 

m 

Spectral  slit  width  due  to  Rayleign  diffraction. 
Electrical  resistivity  (in  ohm-cm) . 

Electrical  conductivity  (in  cgs  esu  unit). 
Dimensionless  conductivity  defined  by  Equation (41) 
Direct  current  conductivity. 

Efficiency  of  the  sphere. 

Relaxation  time. 

Permeability. 

Frequency. 

Dimensionless  frequency  defined  by  Equation  (42) . 
Angular  frequency. 

Plasma  frequency. 

Defined  by  Equation  (56). 
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ABSTRACT 

Hsia,  Jack  Jinn-Goe,  Ph.  D.,  Purdue  University,  August 
1968,  Experimental  Investigation  on  the  Reflectivity  of  Al, 
Cu,  Zn,  and  Ni  and  Its  Comparison  with  Theory.  Joint  Major 
Professors:  Yeram  S.  Touloukian  and  P.  P.  Dewitt. 

Thi-s  study  is  to  generate  accurate  experimental  results 
on  reflectivity  of  metals  and  to  compare  them  with  the  pre¬ 
dictions  made  by  theoretical  models  and  also  to  predict 
values  in  longer  wavelength  range  not  covered  by  this  experi¬ 
ment. 

The  theories  of  the  reflectivity  of  metals  are  first 
reviewed  with  respect  to  macroscopic  and  microscopic  points 
of  view.  Apparatus  was  built  for  measuring  normal  spectral 
reflectivity  by  using  a  multiple  reflection  technique.  Meas¬ 
urements  of  normal  spectral  reflectivity  were  made  on  the 
bulk  surface  of  the  pure  metallic  elements,  copper,  zinc, 
and  nickel  and  high  vacuum  evaporated  surface  of  pure  alumi¬ 
num.  The  measurements  were  made  from  2  to  13/-J  in  air,  at 
room  temperature.  t  '  l-;  (  )( 

The  absorptivity  predicted  by  the  theory  of  anomalous 
skin  effect  assuming  diffuse  electronic  reflection,  one 
electron  per  atom,  and  bulk  conductivity,  is  found  to  have 
lower  values  than  those  obtained  in  the  present  experiment. 

On  the  other  hand,  the  prediction  which  is  made  by  the 


xiii 


simple  equation  derived  from  quantum  theory  is  found  to  fit 
well  with  the  experimental  results. 

Simple  Drude  Theory  and  Bagen-Rubens  relation  are  used 
to  predict  absorptivity  values  in  the  longer  wavelength 
range  not  measured  in  this  experiment.  Most  probable  ab¬ 
sorptivity  values  of  copper,  zinc,  nickel,  and  aluminum 

/ 

from  two  to  one  hundred  microns  are  obtained  based  on 
measurement  and  theory. 


I.  INTRODUCTION 


Attainment  of  better  understanding  of  the  fundamental 
processes  involved  in  the  interaction  of  electromagnetic 
radiation  is  intrinsically  desirable;  further,  it  can  con¬ 
tribute  to  achieving  the  following  four  immediate  objectives. 
First,  a  detailed  knowledge  of  the  mechanism  of  the  inter- 
.  action  of  electromagnetic  fields  with  matter  allows  determin¬ 
ation  of  the  detailed  internal  structure  of  materials  and 
thereby  a  more  complete  knowledge  of  their  physical  proper¬ 
ties.  Second,  if  the  thermal  radiation  properties  can  be- 
related  to  other  more  easily  evaluated  properties  of  mater¬ 
ials,  the  measurement  will  be  simplified.  Third,  if  meaning¬ 
ful  and  practical  equations  can  be  found  or  derived  to  pre¬ 
dict  thermal  radiation  property  data,  obtained  at  a  relative¬ 
ly  few  temperatures  and  wavelengths,  over  much  wider  temper¬ 
ature  and  wavelength  ranges,  the  advantage  is  that  laborious 
and  expensive  measurement  programs  can  be  substantially  re¬ 
duced.  Fourth,  such  a  knowledge  can  help  in  the  development 
and  evaluation  of  equipment  and  techniques  for  measuring  the 
thermal  radiation  properties  of  materials. 

The  first  logical  step  is  to  investigate  the  ideal 
pure  metallic  elements.  Ideal  materials  do  not  need 
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characterization  in  contrast  to  materials  being  used  in  prac¬ 
tical  engineering  environments.  It  is  extremely  difficult  to 
identify  the  surface  characteristics.  Although  ideal  materials 
are  not  of  prime  interest  for  technological  applications,  it 
must  be  appreciated  that  before  real  metallics  or  alloys  can 
be  rationally  studied  and  characterized,  it  is  essential  to 
thoroughly  understand  the  ideal  metals.  This  study  also 
complements  present  activities  at  TPRC  and  in  particular  will 
provide  critical  information  in  support  of  an  analysis  type 
program  to  establish  standard  reference  data  for  the  metallic 
elements . 

Classical  optics  of  metals  has  been  well  reviewed  (12)*, 
and  the  interrelations  of  properties  summarized.  A  signifi¬ 
cant  feature  of  solid-state  physics  during  the  last  two  de¬ 
cades  has  been  the  development  of  the  theory  of  the  anomalous 
skin  effect  in  metals.  This  discovery  has  led  to  a  revolution 
in  ideas  about  the  physical  properties  of  metals,  including 
their  optical  properties  at  low  temperatures.  Quantum  theory 
of  optical  phenomena  in  metals  reveals  much  valuable  informa¬ 
tion  about  their  internal  structure. 

The  range  from  one  or  two  microns  to  longer  wavele  -.gths 
is  the  controlling  region  for  the  heat  transfer  parameters. 

At  450°C  only  1%  of  the  radiant  energy  falls  below  two  mi¬ 
crons  and  at  2000°C  only  about  10%.  The  total  radiation 


* 

Numbers  in  parenthesis  refer  to  the  "List  of  References." 
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properties  are  of  prime  interest  to  the  engineers.  These 
can  be  easily  obtained  by  simply  integrating  the  spectral 
properties  over  the  Wavelengths.  The  modern  technical  trend 
is  toward  studies  at  low  temperatures  and  at  low  temperature 
the  long  wavelength  radiation  properties  are  of  more  practi¬ 
cal  usage.  As  the  first  experimental  phase  of  a  long  range 
program  at  the  Thermophysical  Properties  Research  Center  on 
the  determination  of  thermal  radiation  properties,  an  appar¬ 
atus  capable  of  yielding  high  accuracy  results  has  been  ex¬ 
tensively  evaluated  and  explored.  High  accuracy  of  the 
apparatus  is  required  for  obtaining  precise  data  for  the 
assessment  of  both  the  theoretical  predictions  at  room  tem¬ 
perature  in  the  present  study  and  the  comparison  of  results 
at  different  temperatures  for  other  phases  of  the  program  in 
the  future.  In  the  present  study,  a  multiple  reflection 
technique  is  used.  It  lends  itself  to  high  accuracy  experi¬ 
mental  measurements,  with  the  instrumentation  available, 
normal  spectral  properties  in  the  wavelength  range  of  two 
to  13  microns  are  chosen  for  this  study.  This  is  the  wave¬ 
length  range  of  normal  dispersion  where  simple  theoretical 
equations  are  available. 


II.  THEORIES  ON  REFLECTIVITY  OF  METALS 


The  propagation  of  electromagnetic  waves  within  a 
homogeneous,  isotropic  medium  having  a  conductivity  charac¬ 
teristic  of  a  metallic  conductor  can  be  treated  on  the  basis 
of  the  phenomenological  theory  of  classical  electrodynamics 
(42). 

The  Maxwell's  field  equations  are  expressed  as 


Cur-t  H  =  ~  +  &I2.  E 

cot  c 

Curt  E  =  -  |5 

c  ot 

div  H  =  0 
div  E  =  0 


1 

2 

3 

4 


where  H  is  the  magnetic  field,  E  is  the  electric  field,  eis 


the  dielectric  constant,  a  is  the  electrical  conductivity, 
and  pi  is  the  permeability  which  is  very  close  to  unity  for 
paramagnetic  and  diamagnetic  materials  (42). 

Upon  solving  the  field  equations,  Equations  1  to  4,  the 


electrical  field  can  be  exoressed  as 


E  =  Eq  exp  {i  [■£'  (e  +  ±Z2±)%  - 


where  u)  is  the  angular  frequency  and  z  is  the  penetration 
distance  of  the  wave. 
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The  complex  index  of  refraction  n*  is  defined  as: 


n*  =  (n  -  ik)  -  (E  + 


where/n  is  the  index  of  refraction,  and  k  is  the  index  of 
absorption.  Thus  Equation  5  can  be  written  as 


E  =  EQ  exp  {iW  -  t)}  exp  (-  -^) 


the  phase  velocity  is  reduced  to  c/n, (from  the  first  exponen¬ 
tial  term)  and  the  wave  is  damped  (from  the  second  exponen¬ 
tial  term)  as  it  propagates,  by  a  fraction  exp  (-  2irk)  per 
wavelength. 

Equation  6  yields  the  important  relations  between  the 
electrical  properties  (c,cr)  and  optical  constants  (n,k): 

C  =  n2  -  k2  8 

<7  =  nkco/27T  9 

The  reflectivity,  R,  of  a  metal  is  defined  as  the  ratio 
of  the  radiant  flux  of  the  reflected  and  incident  waves.  The 
radiant  flux  of  the  wave  is  proportional  to  the  square  of  the 
modulus  of  the  electric  field  amplitude.  The  Fresnel  relation 
derived  from  Maxwell  equations  satisfying  continuity  condi¬ 
tions  of  the  electromagnetic  fields  gives  the  well-known 
result  that  for  normal  incidence  (21) 


n*  -  1 
n*  +  1 


6 


or 


R  = 


(n-1) 2  +  k2 
(n+1) 2  +  k2 


11 


The  normal  absorptivity  of  an  opaque  substance  is  related  to 
the  normal  reflectivity  by  the  relation 

R  =  1  -  A  12 


and  by  Equation  11,  the  absorptivity  is  related  to  the 
optical  constants 

A  =  - - - ~  13 

(n+1) A  +  k^ 

For  long  wavelengths  (approximately  X>100  microns) 
(or/y)»e,  where  y  is  oz/2ti .  Equations  8  and  9  yields 

P 

n  =  k  (cr/y)  2  14 

Neglecting  the  higher  order  terms,  -  the  reflectivity  of  a 
metal  from  Equation  11  simplified  to 

R  *  1  -  “  15 

which  leads  to  the  following  well-known'  Hagen-Rubens  equation 
(neglecting  the  higher  order  terms) 


p 

R  =  1-2  (y/cr)  2  16 

where  cris  in  cgs  esu  units.  Equation  16  will  be  in  the  form 
of 

R  «  1  -  36,5  (p/X)*2  17 

if  p  is  in  ohm-cm,  and  X  in  microns. 

As  a  result  of  the  electrical  field  within  the  metal,  a 
current  flow  occurs.  The  current  per  unit  area  is  called 
current  density  and  the  total  current  density  J  can  be  defined 
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(51)  as 

J  =  ~  cur-t  h  16 

Upon  inserting  Equation  1  into  Equation  18  and  replacing  the 
dielectric  constant  e  with  1+4770!  where  a  is  the  polarizabil¬ 
ity,  Equation  18  becomes 


J  - 


1_ 

477 


9E 

St 


+  a 


5e 

<it 


+  CTE 


19 


where  a  and  cr  are  real.  From  Equation  19 ,  it  follows  that 
the  total  current  density  consists  of  three  components; 
namely#  the  pure  displacement  current  (first  term),  the 
polarization  current,  and  the  conduction  current.  The  con¬ 
duction  current  constitutes  the  motion  of  free  charges;  the 
polarization  current  is  associated  with  bound  charges.  Since 
the  pure  displacement  current  is  not  accompanied  by  motion 
of  charges,  the  current  density  may  thus  be  defined  as 

J  =  ol  ||  +  oe  20 


or 


j  =  (or  +  iw  a)E 


21 


It  should  also  be  noted  that  the  phase  difference  between 
the  polarization  current  and  E  is  90°,  so  that  this  current 
does  not  dissipate  the  energy  of  the  field;  the  conduction 
current  is  always  in  phase  with  the  electric  field  and 
therefore  continuously  dissipates  the  energy  of  the  field. 
Hence,  the  absorption  of  electromagnetic  waves  by  a  medium 
is  proportional  with  cr  as  commonly  expressed  by  Ohm's  Law 


J  =  CfE 
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In  the  classical  theory  it  is  assumed  that  the  current 
density  at  a  given  point  is  unambiguously  defined  when  the 
field  strength  at  a  point  is  given.  A  further  assumption 
is  made  that  the  electric  field  may  be  looked  upon  as  uni¬ 
form  when  the  current  density  is  calculated.  This  is  nqt 
valid  when  the  "anomalous  skin  effect"  is  considered  as  will 
be  discussed  later. 


B.  Microscopic  Theory 

1.  Classical  Considerations  -  Simple  Free 
Electron  Model 

The  classical  theory  of  reflectivity  or  absorptivity  of 
metals  based  on  the  free  electron  model  was  developed  by 
Drude  and  Zener* . ■ 

Drude  was  the  first  to  suggest  that  the  fundamental 
physical  properties  of  metals  can  be  explained  by  assuming 
the  existence  of  a  free  electron  gas  moving  between  the  ions 
which  form  the  crystal  lattice  of  the  metal.  These  electrons 
are  in  thermal  equilibrium  with  the  ions.  When  a  constant 
electric  field  is  applied  to  the  metal,  the  electrons  are 
accelerated  in  the  direction  of  the  field.  Therefore,  the 
random  motion  of  the  electrons  is  augmented  by  the  directed 
acceleration,  which  is  responsible  for  the  appearance  of  the 
electric  current.  If  there  were  no  collisions  between  the 
electrons  and  the  lattice,  the  energy  communicated  to  the 
electrons  by  the  field  would  increase  indefinitely. 


See  for  example  in  references  12  and  51. 


resulting  in  an  infinite  mean  free  path,  the  mean  distance 
travelled  by  electrons  between  collisions.  This  would  lead 
to  zero  electrical  resistivity.  The  fact  that  metals  have 
a  finite  resistivity  suggests  that  collisions  between  elec¬ 
trons  and  lattice  ions  do  occur. 

When  the  damping  due  to  the  resistance  of  the  metal  is 
considered,  the  equation  of  motion  for  a  free  electron  is 


mr  +  myr  =  eE  23 

where  y  is  the  damping  coefficient,  r  is  the  Electron  posi¬ 
tion  coordinate  within  the  metal;  E  =  EQe  is  the  applied 
electromagnetic  field;  and  e  and  m  are  the  charge  and  mass 
of  an  electron. 

Upon  solving  the  equation  of  motion  (Equation  23)  for 
the  velocity  r,  the  current  density  J  is  obtained  by  using 
the  relation  J  -  Ner.  The  electrical  properties  (ne )  are  ob 
tained  from  the  relation  between  the  current  density  and 
electric  field  (Equation  21).  Equations  8  and  9  for  the 
relations  between  the  electrical  properties  (cr,e)  and  opti¬ 
cal  constants  (n,k)  then  give  the  Drude-Zener  formulas: 


2  ,  2  ,  477Ne‘ 

n  -k  =  1  -  — — — 


nk  = 


277  Ne 
mco 


2  2 

co  +y 


2  ,.2 

co  +y 


24 

25 


where  N  is  the  number  of  electrons  per  unit  volume .  The  co¬ 
efficient  y  is  related  to  tne  dc  electrical  conductivity  by 


the  relation 


10 


v  -  M 
'  ~  mcr 


2 
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The  normal  spectral  reflectivity  or  absorptivity  can  be  ex¬ 
pressed  in  terms  of  e,  m,  y  by  means  of  Eauations  11,  13, 
24,  and  25. 

The  results  of  the  Drude-Zener  theory  when  compared  with 
experimental  measurements  show  fair  agreement  for  the  liquid 
metals  Hg  and  Ga  (31),  and  ultra-high  vacuum  deposited  silver 
gold,  and  aluminum  (4).  In  general,  the  theory  fails  in  the 
higher  frequency  region  of  the  spectrum  and  also  at  low  tem¬ 
peratures  . 


2.  Multi-free/bound  Electron  Theory 
Roberts  has  suggested  (44,45,46)  and  the  quantum  theory 
has  also  exposed  the  idea,  that  u  real  metal  consists  of 
several  types  of  free  and  bound  electrons  existing  in  a  metal 
and  both  types  of  electrons  exist  in  different  energy  states 
(12).  As  such,  they  react  in  different  degrees  to  the  in¬ 
fluence  of  an  electric  field  induced  by  the  incident  radi¬ 
ation.  The  free  electrons  of  a  particular  type  are  distin¬ 
guished  by  the  damping  coefficient  as  previously  suggested. 
The  bound  electrons  influenced  by  the  induced  electric  field 
can  be  assumed  to  be  characterized  by  an  elastic  restoring 
force  and  a  viscous  damping  force.  Equations  of  optical 
properties  and  radiative  properties  can  be  derived  by  start¬ 
ing  with  the  equations  of  motion  for  different  types  of 
electrons  using  the  procedure  discussed  in  the  previous 
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section.  The  resulting  equations  are  rather  complicated  and 
involve  parameters  mentioned  above  as  vf^ll  ?s  the  number 
density  of  the  different  types  of  electrons. 

Roberts  has  tested  this  model  for  several  metals  and 
evaluated  the  constants  required  from  measurements  of  the 
optical  constants.  In  each  ca~e  of  the  metal  observed,  it 
was  possible  to  determine  values  of  a  finite  number  of  the 
constants  which  gave  an  acceptable  comparison  between  data 
and  theory.  Attempts  were  made  by  Richmond,  et  al. (41)  to 
determine  the  values  of  the  parameters  which  would  fit  the 
normal  spectral  reflectivity  spectrum  of  rhodium  by  using 
the  multi-free/bound  electron  model.  The  comparison  of 
the  theory  and  data  was  very  favorable.  However,  the  opti¬ 
cal  constants  computed  from  the  parameters  in  the  visible 
portion  of  the  spectrum  (where  the  computed  curve  shows  good 
agreement  with  experimental  data)  did  not  agree  with  meas¬ 
ured  optical  constants  for  the  same  wavelengths.  Hence,  it 
is  evident  that  this  model  suffers  from  its  cumbersome  form 
and  large  number  of  parameters,  giving  an  appearance  of 
mathematical  expediency  rather  than  physical  reality.  Addi¬ 
tional  tests  of  this  model  require  optical  constants  data  or 
complete  reflectivity  spectra.  For  testing  this  model  in 
the  infrared  range  the  parameters  for  the  bound  electrons 
can  be  dropped  because  they  are  not  important  in  this  wave¬ 
length  range. 
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3.  The  Anomalous  Skin  Effect 
a.  Basic  Concept.  As  the  temperature  is  lowered,  the 
depth  of  penetration*,  6 , of  the  high  frequency  field  into  the 
metal  decreases  as  a  result  of  the  increase  in  electrical 
conductivity;  the  mean  free  path  increases  and  may  become 
comparable  with  or  much  greater  than  the  depth  of  penetra¬ 
tion.  In  this  condition  the  electrons  will  succeed  in 
traversing  regions  with  different  field  strength  in  one  mean 
free  path,  and  the  additional  velocities  which  they  will  ac¬ 
quire  as  a  result  of  this  process  will  depend  on  the  field 
strength  along  the  entire  path.  This  means  that  Ohirrs  Law 
J  =  <TE  in  which  a  is  constant  for  all  parts  of  the  metal, 
must  be  replaced  by  a  more  general  equation  in  which  J  should 
be  a  function  of  field  strength  and  position.  This  is  the 
anomalous  skin  effect  in  which  the  expression  of  the  electric 
field  strength  is  given  by  a  very  complicated  form,  derived 
from  the  electron  distribution  function  and  Maxwell  field 
equations,  and  is  not  exponential  in  form  as  in  the  case  of 
classical  theory  (refer  to  Equation  7).  As  long  as  the  pro¬ 
pagation  of  the  wave  is  no  longer  exponential  in  nature,  the 
classical  representation  involving  the  complex  refractive 

index,  defined  from  the  exponential  form  of  the  electric 
★ 

Referring  to  Equation  7,  the  depth  of  penetration,  defined 
as  the  distance  traversed  by  the  wave  within  the  material 
where  the  wave  amplitude  is  reduced  by  a  factor  of  1/e,  can^ 
be  expressed  as  6  =  A/27rj<.  At  long  wavelengths,  k  =  {oq/v)  , 
therefore  6  =  c/ (2ituj  <7q) This  is  not  to  be  confused  with 
the  skin  depth,  6S  =  (mc^^TrNe2)^,  which  is  the  distance  of 
one  wavelength  corresponding  to  the  plasma  frequency. 
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field,  is  no  longer  physically  significant.  The  more  general 
term  "surface  impedance"  2  is  introduced  and  defined  as  the 
ratio  of  the  electric  field  at  the  surface  of  the  metal  to 
the  total  current  per  unit  area.  By  using  Maxwell’s  field 
equations,  the  surface  impedance  Z  can  be  expressed  as 


(E) 


2=0 


'  aEN 
dz J z=0 


27 


or 


2  -  4»c  <E)z=0/(H)2=0 

At  normal  incidence  the  absorptivity  in  terms  of  the 
impedance  concept  is  given  by 


28 


R  =  1  -  A  = 


4m /c  -  Zl 


|4m/c  +  zl 

where,  4m /c  and  Z  are  the  impedance  of  free  space  and  metal 

surface  respectively.  Equation  29  can  also  be  written  as 

2 


29 


R  =  1  -  A  = 


4 m  1  , 

c  Z 


4m  JL 
c  Z 


+  1 


30 


Comparison  between  Equations  30  and  10  reveals  the  re¬ 
lationship  between  the  surface  impedance  and  complex  index 
of  refraction  as 


31 


Furthermore,  when  the  surface  impedance  Z  is  expressed 
as  r  +  ix,  where  r  is  the  surface  resistance  and  x  is 
the  surface  reactance,  expanding  Equation  29  in  series  yields 


*-?{*-£(?)  —Mdr) 


x2 
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A  =  ££. 
A  IT 


The  second  and  higher  powers  of  x  may  be  neglected.  This 

is  tantamount  to  neglecting  the  displacement  current,  the 

$  .... 
photoelectric  current  ,  and  the  atomic  polarizability  since 

these  currents  are  not  dissipating  the  energy  of  the  electric 

field  just  as  the  condition  of  the  surface  reactance  x.  The 

term  r  can  also  be  neglected  if  the  absorbing  power  is 

small,  a  fact  which  is  valid  for  good  conductors.  Thus  the 

following  final  expression  for  the  absorptivity  is  obtained. 

33 

The  microscopic  phenomena  of  metals  upon  receiving  an 
electromagnetic  wave  may  be  discussed  as  follows.  The  con¬ 
duction  electrons  in  the  metal  interact  with  the  lattice  and 
with  each  other.  An  electromagnetic  wave  is  incident  on  the 
surface  of  the  metal.  The  electron  energy  distribution 
(Equilibrium  Fermi  distribution)  is  modified  under  the  action 
of  the  electromagnetic  field  and  as  a  result  of  collisions 
between  electrons  and  between  electrons  and  the  lattice.  If 
the  change  in  the  equilibrium  function  is  small,  the  pertur¬ 
bation  theory  can  be  used  to  obtain  the  electron  distribution 
function,  f,  then  the  current  density  can  be  calculated  by 


34 


using  the  following  Equation  (37). 

3 

J(z)  =  -2e  (a)  JJJvxfdvxdvydv2 

The  electric  field  strength  can  be  obtained  by  using  the 
following  equation  derived  from  the  Maxwell's  field  equations. 


*The  photoelectric  current  is  due  to  the  flow  of  electrons 
emitted  fcom  the  surface  as  a  result  of  the  interaction  be¬ 
tween  the  conduction  electrons  and  the  electromagnetic  field. 
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+  %E  =  -^4^  J  *  35 

dz^  c  c 

Once  the  expression  of  the  electric  field  strength  is  derived, 
the  surface  impedance  Z  can  be  obtained  by  using  Equation 


b.  Normal  Skin  Effect.  The  theory  of  normal  skin  effect 
(which  is  identical  to  the  classical  theory  discussed  in  the 
previous  section  as  will  be  seen  in  the  following  develop¬ 
ment)  differs  from  the  theory  of  the  anomalous  skin  effect. 
In  the  former  theory  the  change  of  the  electron  distribution 
function  in  the  penetration  direction  is  neglected,  a  fact 
which  is  valid  when  the  mean  free  path  is  much  smaller  than 
the  depth  of  penetration.  The  electric  field  strength  thus 
derived  is  exponential  in  form  and  the  surface  impedance  is 
given  as 


2  [Sit  _  .1 -  'j 
2  12  2  ..  J 


C  uc  C  (1+ico  T; 

where  r  is  the  relaxation  time  which  is  equal  to  y  Using 
the  relation  between  the  surface  impedance  and  complex  index  of 
refraction  as  stated  in  Equation  31,  Equation  36  has  the 
same  form  as  the  Drude-Zener  formulas  (i.e..  Equations  24 
and  25  when  Ne  /m  is  replaced  by  yo  as  stated  in  Equation  26) . 

At  sufficiently  long  wavelengths  (co  T«l)  the  free  elec¬ 
trons  dominate,  thus  the  displacement  can  be  neglected  since 
it  is  not  associated  with  charges.  Equation  36  becomes 

Z  =  r  +  ix  =  /~^  (1+i)  37 
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From  Equation  37  the  surface  resistance  is  expressed  as 


38 


c.  Anomalous  Skin  Effect.  When  the  mean  free  path  is  com¬ 
parable  with  or  much  greater  than  the  depth  of  penetration, 
the  change  of  the  electron  distribution  function  in  the  pene¬ 
tration  direction  can  no  longer  be  neglected.  With  the 
above  consideration  Reuter  and  Sondheimer*  developed  the 
quantitative  theory  of  the  anomalous  skin  effect.  Dingle 
(14)  simplified  their  expressions  in  the  visible  and  infra¬ 
red  region  (0.1<X<30  microns)  and  gave  the  expressions  of 
absorptivity  for  the  following  two  special  cases:  (1)  per¬ 
fectly  specular  reflection  of  electrons  at  the  inner  surface 
of  the  metal  (p  =  1)  and  (2)  perfectly  diffuse  reflection 
(p  =  0) ,  namely. 


and 
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40 


where  v^  is  the  Fermi  velocity  and  a  and  v  are  the 
dimensionless  conductivity  and  frequency,  respectively,  given 
by 


See  for  example  in  reference  51. 
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rvf/ 

a  =  arL-^ 


-  =  •'[I  e(rt)  ] 


V-  vf  -Ne*.-  - 

Holstein  (31),  'sing  energy  balance  considerations, 
derived  the  absorptivity  expression  for  the  case  when  the 
mean  free  path  approaches  infinity,  i.e.,  at  very  low 
temperature 

2lTHe2v.3  3v. 

A  =  P  - 2 -T~  +  I1  '  P>  ST  43 

m  w  c 

where  p  is  the  fraction  of  electrons  which  are  reflected 
specularly.  Thus,  for  the  case  of  diffuse  reflection  of 
electrons 

3  vf 

A0  =  4  ~  44 

And  for  the  case  of  specular  reflection 
2ffNe2vf3 

A1  =  2  3  43 

mo:  c 

Expressed  in  terms  of  the  dimensionless  frequency  as  stated 
in  Equation  42,  Aj  can  be  written  as 

2  vf  — ~2 

A1  "  3  —  ^  46 

Equation  46  can  also  be  obtained  from  Equation  40  when 
allowing  a  to  approach  w,  as  is  the  case  for  very  low 
temperatures. 


4.  Quantum  Theory 


From  Equations  9  and  25  of  the  classical  theory,  the 
conductivity  .can  he  expressed  as 


Ne _ y 

°  ~  ra  2  2 

co  +  y 


u  4tt  2  ,  2 

co  +  y 

2  h 

where  =  (47rNe  /m)  is  the  frequency  characteristic  of  the 
electron  plasma  oscillation.  The  plasma  frequency  is  gener¬ 
ally  outside  the  high-frequency  limit  of  the  free-electron 
dispersion  spectrum.  For  example,  the  plasma  frequencies  of 
silver,  gold,  copper,  and  aluminum  are  0.31,  0.50,  0.575, 
and  0.8  micron,  respectively  (34,37). 

Holstein  (31)  used  the  perturbation  theory  to  obtain  the 
damping  coefficient  y.  With  y  known  as  a  function  of  wave¬ 
length  and  temperature,  the  conductivity  a  can  be  calculated 
by  using  Equation  48.  The  dielectric  constant  e  should  be 
obtained  from  the  Kraners-Kronig  relation 


—  f 

ir  J 


‘iTJ  7T  J  n  /  i  2  2. 

0  (co  -co  ) 

since  upon  relating  Equations  24  and  8,  the  expression  for 

c  obtained  in  the  form 
2 

1-e  1 

47 T  477  2  .  2 

co  +y 


which  is  inconsistent  with  the  result  when  a  in  the  expres¬ 
sion  of  Equation  47  is  inserted  in  Equation  49.  The  reason 
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for  such  inconsistency  is  due  to  the  fact  that  in  quantum 
theory,  the  damping  coefficient  y  in  the  expression  of  a 
(Equation  47)  is  a  function  of  not  only  temperature  but  also 
wavelength.  However,  in  the  classical  theory,  y  is  a  func¬ 
tion  of  temperature  only.  Therefore,  Equations  49  and  50  can 
be  used  interchangeably  in  the  qlassical  theory.  Once  the 
electrical  properties  (a,e)  are  known,  the  optical  constants 
(n,k)  can  be  derived  from  Equation  8  and  9.  Furthermore, 

the  absorptivity  is  obtainable  by  using  Equation  13  which 

\ 

relates  the  absorptivity  and  the  optical  constants.  Some  of 
the  considerations  that  enter  into  the  calculations  are 
a.  free  electron  theory,  b.  spherical  Fermi  surface,  c. 
Fermi-Dirac  Statistics  of  the  electron  distribution,  d.  Pauli 
exclusion  Principle,  and  e.  the  classical  treatment  of 
electromagnetic  field. 

It  is  well-known  that  the  power  expenditure*  W  for  a 
free  electron  system  is  related  to  the  conductivity  cr  (a;)  by 
the  relation: 


NW  =  ~cr  (W)  51 

where  W  is  the  power  expenditure  due  to  one  electron  per 

2 

unit  time  and  E  /2  is  the  ener'/y  density  per  unit  volume 
of  the  electromagnetic  field. 


W  is  the  time  average  of  -W  =  (dS  /hz)  where  W  is  the 
heat  produced  per  unit  time  and  pe§  unit  volume  within 
conducting  medium,  z  is  the  distance  of  penetration  and 
is  the  Poynting  vector  which  is  defined  as 

s  =  ^  (B  X  a) . 


Joule 
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The  power  expenditure  W  is  derived  from  the  transition 
probability  which  is  obtained  by  using  the  perturbation  theory 
on  solving  the  Schrodinger  equation  for  a  conduction  electron 
interacting  with  the  perturbing  electromagnetic  field  and  the 
lattice  vibration. 

The  damping  coefficient  yQ^  due  to  electron-phonon 
interaction  is  given  as: 

ye  p(V'a)  =  yepla)hep(ftf0t)  52 

Here 


where  It  is  Planck's  constant/2ir,  K  is  Boltzmann  constant, 
©P  is  the  Debye  temperature,  and  ^ep(M»a)  *-s  t^ie  quantum 
correction  factor  due  to  electron-phonon  interaction. 

The  damping  coefficient  including  the  effect  due  to 
electron-electron  interaction  is  given  as 

:  y(n*a)  =  yep(a)bep^,0!)  +  Vee^'a*  53 

where  yee(l*<0')  is  the  damping  coefficient  due  to  electron- 
electron  interaction. 

2  2  2 

In  the  spectral  region  y  «oo  «u^  ,  Holstein  obtained 
the  following  equations  for  optical  constants  (n,k): 


n 


=  2~2  ^ep(0!)bep 


(a)  + 


JsL 


(a)‘ 


In  this  region  the  quantum  correction  factor  bep  is 


54 

55 
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independent  of  wavelength  and 

“  (I?)2  £<“> 

Thus,  the  bulk  absorptivity  An  can  be  expressed  as 

13 


Notice  from  Equation  54  that  n  is  proportional  to 

2 

(A  +  constant)  and  k  is  proportional  to  A  for  a  fixed 
temperature.  Thus,  the  plottings  of  n  vs  A  and  k  vs  A 
will  give  or,  0  and  y®  (a)  since  b  (a)  is  known  from 

the  theory.  The  quantity  (°0  can  be  calculated  by  Equa 

tion  56.  With  all  these  parameters  known,  the  bulk  absorp¬ 
tivity  can  be  easily  obtained  from  Equation  57.  Direct 
current  conductivity  can  be  calculated  by  using  the  equation 
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III.  METHODS  OF  NORMAL  SPECTRAL  REFLECTIVITY 

MEASUREMENTS 


A.  Calorimetric  Method 

The  basic  principle  is  to  find  the  absorptivity  by  a 
comparison  of  the  energy  absorbed  with  the  energy  of  the 
normally  incident  radiation  in  the  absolute  method  or  to 
find  the  absorptivity  relative  to  a  black  surface  in  the 
comparative  method. 

1.  Absolute  Measurement 

Biondi  (7)  used  the  calorimetric  method  to  determine 
the  absorptivities  of  copper  and  silver  at  3.4  to  4.2°K 
over  the  wavelength  range  from  0.3  to  3.3  microns..  The 
energy  of  the  incident  radiation  absorbed  on  the  surface 
of  a  metal  specimen  is  determined  by  first  measuring  the 
rise  in  temperature  of  the  specimen  and  then  by  turning 
off  the  incident  radiation  source  and  power  being  applied 
to  the  heater  in  the  specimen  stage  (see  Figure  1)  to 
reproduce  the  temperature  rise.  The  energy  applied  to 
the  specimen  heater  is  equivalent  to  the  incident  energy 
absorbed  by  the  specinu3  .  The  total  incident  energy,  the 
energy  absorbed  by  the  specimen  plus  the  energy  reflected  by 
the  specimen  and  absorbed  by  the  black  absorber,  is 


Specimen 

Specimen  Stage 

(Thermometer 
&  heater) 

Copper  base 
&  thermometer 


Figure  1.  Schematic  of  Biondi's  Apparatus 
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equivalent  to  the  energy  required  to  put  into  the  heaters  in 
both  the  specimen  stage  and  the  absorber  stage  to  reproduce 
the  temperature  rise.  The  temperature  was  determined  by  a 
carbon  resistance  thermometer.  The  precision  was  claimed 
around  +2%. 

2.  Comparative  Measurement 

Brandenberg,  Clausen,  and  McKeown  (8)  measured  the  ab¬ 
sorptivity  of  evaporated  aluminum  in  the  wavelength  region 
from  0.45  to  2.0  microns  by  using  a  specimen  which  consists 
of  a  temperature-sensitive  Y-cut  quartz  crystal  plated  on  one 
side  with  the  metal  film  to  be  studied  and  the  other  side 
with  an  opaque  layer  of  benzene  smoke,  a  highly  absorbing 
coating.  The  small  temperature  rise  of  the  plated  crystal 
due  to  the  energy  absorbed  from  the  incident  beam  of  radi¬ 
ation  is  measured  by  monitoring  the  frequency  change  of  the 
quartz  crystal  when  it  is  used  as  an  active  element  and 
driven  in  an  oscillator.  The  absorptivity,  relative  to  the 
black  surface,  is  determined  from  frequency  rates  at  a  par¬ 
ticular  temperature  as  the  crystal  goes  through  a  heating 
and  cooling  cycle.  This  method  is  claimed  to  give  a  possible 
precision  in  absorptance  of  jd.,4%.  The  schematic  of  the 
apparatus  is  shown  in  Figure  2. 

The  crystal  holder  is  mounted  rigidly  to  a  metal  hemis¬ 
phere  such  that  the  center  of  the  crystal  coincides  with  the 
origin,  and  the  black  face  is  subtended  by  the  hemisphere. 
This  increases  the  blackness  of  the  benzene  smoke  by 


Monochromator 


Figure  2.  Schematic  of  the  Apparatus  Used 
by  Brandenberg,  et  al. 
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refocusing  reflected  energy  back  onto  the  benzene  surface. 

The  crystal  can  be  turned  through  180°  so  that  the  beam  is 
either  incident  on  the  sample  electrode  or  incident  through 
the  hemisphere  on  the  benzene  smoke  coating . 

B.  Single  Reflection  Method 

Radiation  energy  received  from  a  single  reflection  from 
the  specimen  is  compared  with  the  direct  incident  energy. 
Twidle's  apparatus  (53)  illustrates  this  method  (see  Figure 
3).  The  mirror  (specimen)  under  investigation,  M,  could  be 
,Amng  into  the  beam  so  that,  after  reflection,  the  image  of 
the  slit  S  is  formed  on  the  thermopile  T2»  Thus,  this  appar¬ 
atus  requires  different  detectors  to  measure  the  incident  and 
reflected  radiant  flux.  This  will  give  rise  to  large  errors 
in  the  measurement  results. 

C.  Multiple  Reflection  Method  (Strong's  Method) 

Multiple  reflection  techniques  have  been  discussed  by 
Strong  (52)  and  Beattie  (1).  The  basic  principle  of  Strong's 
method  can  be  illustrated  by  the  original  setup  proposed  by 
Strong  (see  Figure  4) .  The  square  of  the  absolute  reflec¬ 
tance  of  the  sample  can  be  obtained  from  the  ratio  of  two 
radiant  flvx  readings  corresponding  to  the  two  conditions 
shown  in  Figure  4.  Because  the  square  of  the  absolute  reflec¬ 
tivity  is  measured,  the  value  of  the  absolute  reflectivity  is 
obtained  with  higher  precision  than  would  be  possible  of 
there  was  only  one  reflection  from  the  specimen. 
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Figure  3.  Schematic 


S  =  Exit  Slit 
Mc  Specimen 
T,,T2b  Thermopiles 

of  Twidle's  Apparatus. 
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Standard  mirror,  r* 


Figure  4.  The  Strong  Technique. 


1.  Gates,  Shaw,  and  Beaumont  (22)  Method 
Gates  et  al.,  used  the  modified  Strong's  method  by  meas¬ 
uring  the  reflectivity  of  a  matched  pair  of  evaporated  speci¬ 
mens,  size  three  inches  by  ten  inches.  The  flux  of  the  radi¬ 
ation  after  being  reflected  from  the  parallel  specimens  2,4, 
6,8,  and  10  times  was  recorded,  and  the  reflectivity  was 
determined  as  a  root  of  the  ratios  of  these  readings.  The 
incident  angles  were  20°  to  60°,  wavelength  range  0.7  to  12 
microns,  and  the  specimens  consisted  of  evaporated  metal 
films.  The  accuracy  of  the  measurement  was  said  to  be  +0.2%. 
(Collimated  radiation  was  used  as  the  source.) 

2.  Fowler  (Far  Infrared)  Method  (20) 

Fowler  also  used  the  modified  Strong's  method  in  the 
wavelength  range  8.5  to  83.5  microns.  The  reflectivity  was 
calculated  from  the  ratio  of  the  observed  radiant  flux  after 
seven  reflections  to  that  observed  after  one  reflection  from 
the  specimens.  The  incident  angle  was  23°,  and  the  matched 
pairs  of  specimens  were  evaporated  gold  films  on  glass  sub¬ 
strates.  The  specimen  sizes  were  7.6  cm  X  3.8  cm  and  4.8  cm 
X  3.8  cm.  The  accuracy  stated  was  +0.2%. 

3.  Bennett  (Visible,  Averaging  Sphere)  Method (6) 

Strong ' s  method  was  used  with  both  specimen  and  standard 
mirrors  included.  An  averaging  sphere  was  used  to  reduce  the 
errors  introduced  by  optics  misalignment  and  spatial,  angular 
and  area  sensitivities  of  the  detectors.  Two  block-type 
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specimen  holders  were  used.  Measurements  were  made  with 
upper  and  lower  blocks  interchanged  to  eliminate  the  differ¬ 
ence  of  the  standard  mirrors  and  also  to  reduce  the  geometri¬ 
cal  defect  of  the  blocks  (e.g.,  being  not  exactly  parallel). 
The  angle  of  incidence  is  7°,  with  two  reflections  on  the 
specimen,  and  the  accuracy  is  +0.1%. 

4.  Bennett  (Visible  and  Infrared)  Method  (6) 

About  1960,  materials  for  the  inner  coating  of  the  aver¬ 
aging  sphere  to  be  used  in  the  near  infrared  range  were  not 
explored.  Thus,  Bennett  developed  an  optical  system  with  a 
unique  compensating  feature  which  prevented  the  image  on  the 
detector  from  changing  size  or  position  because  of  a  slight 
tilt  of  the  specimen.  The  angle  of  incidence  is  5°,  with 
two  reflections  from  the  specimen.  The  specimen  size  is  1.5- 
inch  diameter  disk.  The  accuracy  is  +0.1%. 

5.  Dunn  (infrared.  Averaging  Sphere)  Method  (17) 

Dunn  used  averaging  spheres  to  reduce  the  optical  system¬ 
atic  errors.  He  used  sulfur  as  the  averaging  sphere  coating 
(18)  which  has  a  reflectance  of  approximately  0.95  to  0.75 
for  the  wavelengths  ranging  from  two  to  seven  microns,  re¬ 
spectively.  The  modified  Strong's  method  was  used  to  measure 
the  reflectance  of  matched  pairs  of  specimens  with  size  of 
0.5  inch  X  0.25  inch.  The  optical  path,  when  the  specimens 
are  in  and  out  of  the  light  path,  were  not  identical.  The 
incident  angles  were  0  to  52°,  with  one  reflection  of  each 
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specimen  of  the  matched  pairs.  The  accuracy  was  given  as 
+0.15%. 


D.  General  Comments 

in  the  previous  sections,  various  methods, such  as  the 
calorimetric  method,  single  reflection  method,  and  multiple 
reflection  method,  have  been  discussed. 

Although  the  calorimetric  method  gives  high  accuracy 
in  reflectivity  (for  low  absorptivity  materials),  it  is  dif¬ 
ficult  to  use  the  calorimetric  methods  in  the  infrared  be¬ 
cause  of  the  lo*-.  flux  of  the  monochromator  radiation  avail¬ 
able.  The  comparative  measurement  of  calorimetric  method 
appears  capable  of  detecting  very  small  temperature  changes 
in  the  specimen  only  when  the  specimen  is  an  evaporated  film 
on  Y-cut  quartz  crystal.  Therefore,  bulk  specimens  cannot 
be  measured  by  using  this  method. 

Single  reflection  method  needs  two  detectors  to  receive 
incident  and  reflected  radiation,  and  the  accuracy  of  the 
measurement  is  poor.  On  the  other  hand,  multiple  reflection 
method  generally  yields  more  accurate  results  than  the  single 
reflection  method.  In  the  methods  used  by  Gates  et  al.,  and 
Fowler,  the  radiant  beam  undergoes  more  than  two  reflections 
on  the  specimen  surface,  thus  a  relatively  large  size  speci¬ 
men  is  required  which  will  increase  the  difficulty  in  prepar¬ 
ing  the  bulk  specimen.  Also, in  these  apparatus  because  the 
incident  angle  has  to  be  large,  the  measurement  can  no  longer 
be  considered  as  a  normal  reflectivity  measurement  (6) . 
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Bennett's  apparatus  in  both  the  visible  and  the  infrared 
range  and  Dunn's  apparatus  in  infrared  range  all  employed 
special  features  to  reduce  the  optical  systematic  errors. 


.£V.  PRINCIPLE  OF  MEASUREMENT 


An  apparatus  for  measuring  absolute  normal  spectral 
reflectivity  is  desired.  Some  special  requirements  are: 
the  apparatus  should  be  able  to  measure  reflectivity  in  the 
infrared  range  from  two  to  13  microns,  the  apparatus  should 
give  accurate  results  (jf0.002  reflectivity  unit),  the  inci¬ 
dent  angle  on  the  specimen  should  be  less  than  10°,  the 
apparatus  should  be  able  to  measure  both  thin  film  and  bulk 
specimens,  the  specimen  size  required  should  be  small  (1x2 
inches  or  1-1*4  inch  diameter)  to  simplify  preparation  of  bulk 
specimens,  and  that  optical  systematic  errors  are  kept  as 
small  as  possible. 

With  these  requirements  in  mind  and  from  the  comments 
of  the  previous  section,  the  logical  selection  is  an  appar¬ 
atus  employing  the  multiple  reflection  method,  an  apparatus 
similar  to  that  of  Bennett's  using  two  identical  blocks  as 
specimen  holder  and  employirg  a  flux  averaging  sphere  with 
an  inner  coating  usable  in  the  infrared  range. 

Following  is  the  principle  of  operation  of  the  block 
type  specimen  holder.  The  measurements  are  made  for  four 
orientations  of  the  specimen  (see  Figure  5)  in  order  to  re¬ 
duce  the  effects  of  different  reflectivities  of  standard 
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mirrors  and  to  compensate  for  the  geometrical  effects  of  the 
block  surfaces.  By  taking  the  products  of  the  ratios  meas¬ 
ured  in  four  orientations,  the  absolute  reflectivity  of  the 
specimen  can  be  obtained  even  if  the  reflectivities  of  the 
standard  mirrors  mounted  on  the  two  blocks  are  not  exactly 
equal. 

For  1%  noise  level  of  the  output  signal,  the  uncertainty 
of  the  final  results  of  the  reflectivity  measurements  based 
on  the  radiant  flux  ratio  can  be  estimated  as  follows: 

For  each  of  the  orientations  of  the  blocks,  reflectivity  is 
given  by 


Ri  =  (I"/!')*2 

Following  Kline  and  McClintock  (35) ,  for  uncertainties 
in  single-sample  experiments,  the  uncertainty  of  R. ,  A  R^,  is 


For  Ri  **  1  I "  sb  1 1 


A  R. 


R 


Zk  =  / 1 I  Ail 

LL2  i"  J 


,2 


and  for  R^  «  1, 


Ail  =  Ai 

i'  I 


therefore 
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=  2“^  (4~)  =  0.007 

since  1%  noise  level  has  been  assumed  ^i.e.  =  0.01^). 

The  final  results  of  reflectivity  is  calculated  from 
observations  on  four  different  orientations.  It  is  expressed 
as 

R  =  (R^R^ 

Following  Kline  and  McClintock  (35) ,  the  uncertainty  of 
R,  A  R,  is 


since  it  is  quite  safe  to  assume  that  each  term  in  parenthesis 
gives  the  same  value. 

Thus 


a 

A  R  =  2  *  ~  (r1R2R3R4)4  R2R3R4A  R1 


A  R 
R 


1  A  ?1 

2  '  R, 


A  R 
R 


/  A  \ 

=  0.0035  (^since  — —  =  0.007y 


A  R  =  0.0035  when  R  1 

Thus,  the  uncertainty  of  each  data  point  is  expected  to  be 
about  0.0035  reflectivity  unit  assuming  an  estimated  uncer¬ 
tainty  in  the  output  signal  measurement  of  1%.  in  reality 
this  is  larger  than  the  observed  uncertainty  which  falls  with¬ 
in  +0.002. 
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V.  DESCRIPTION  OF  THE  APPARATUS 


A.  General  Features 


An  apparatus  was  designed  and  constructed  to  measure  the 
normal  spectral  reflectivity  by  using  the  multiple  reflection 
technique  in  the  infrared  region,  2  to  13  microns,  in  air  and 
at  room  temperature.  The  angles  of  incidence  and  reflection 
are  fixed  at  seven  degrees  from  the  normal  to  the  specimen 
surface.  The  apparatus  consists  of  the  conventional  Perkin- 
Elmer  12-C  single  beam  single  pass  infrared  spectrometer, 
transfer  optics,  specimen  holder,  and  recorder.  The  block 
diagram  and  general  arrangement  of  the  system  are  shown  in 
Figures  6  and  7,  respectively. 


B.  Specimen  Holder 

The  specimen  holder  is  composed  of  two  identical  blocks, 
A  and  B,  made  of  stainless  steel  304.  The  overall  dimensions 
for  each  block  is  3.75  X  3.75  X  2.0625  inches.  Two  cut-outs 
are  made  for  the  light  path.  The  upper  and  lower  surfaces 
of  each  block  are  machined  parallel  to  each  other  wichin 
0.0005  inch.  Both  the  inside  and  the  outside  walls  of  the 
blocks  are  sprayed  with  3M  velvet  black  coating  to  minimize 
reflected  stray  energy.  The  blocks  are  connected  to  each 
other  with  two  positioning  pins,  and  either  of  which  may  be 


Figure  7.  Schematic  of  Apparatus  For  Measurement  of  Absolute  Normal 

Spectral  Reflectivity. 
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placed  on  top  of  the  other.  The  two  standard  mirrors 
(Electroplated,  first  surface  rhodium  plane  mirrors,  size 
1/2  x  1  3/8  inches.  Evaporated  Metal  Films  Corporation, 
Ithaca,  New  York)  and  the  specimen  under  study  are  mounted 
by  clip  springs  on  the  surfaces  of  the  blocks.  The  rhodium 
mirrors  are  chosen  because  their  reflectance  values  in  air 
are  quite  stable  (29,43).  The  specimen  holder  and  the  loca¬ 
tions  of  specimen  and  standard  mirrors  are  shown  in  Figures 
8  and  9. 

The  specimen  holder  can  be  manually  slid  on  a  supporting 
base  to  positions  of  SPECIMEN  IN  and  SPECIMEN  OUT.  The  posi¬ 
tions  are  located  by  stops  at  the  two  ends  of  the  supporting 
base  with  adjustable  screws  for  alignment.  The  supporting 
base  is  in  turn  supported  by  three  brass  columns  which  are 
fastened  on  a  leveling  base  which  can  be  adjusted  by  three 
pairs  of  screws  for  optical  alignment. 

C .  The  Transfer  Optics 

Figure  7  and  Table  1  indicate  the  transfer  optics,  their 
dimensions,  and  relative  distances.  All  transfer  optics  are 
first  surface  aluminum  mirrors.  The  exit  slit  is  imaged  on 
the  standard  mirror  on  the  specimen  holder  by  1:2  ratio  and 
re-imaged  by  1:1  ratio  on  the  detector.  The  average  angle 

A 

of  incidence  on  the  surface  is  about  7~.  The  horizontal 
folding  of  the  infrared  beam  is  necessary  because  of  the 
limited  space  of  the  base  plate  of  the  optical  bench  on 
which  all  the  mirrors,  monochromator,  and  the  specimen 
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Table  1.  Transfer  Optics 


Mirror 

Kind 

Sizi  (in. ) 

Focal 

Length (in. ) 

Distance 
(in. ) 

M]_ 

s<*> 

6 (dia . ) 

16 

10  (d) 

m2 

D(b) 

4x4 

21 

M3 

D 

3x3 

10 

M4 

D 

lfcxl 

14<e) 

m5 

D 

1*4X1^ 

15 

M6 

S 

4-V(dia. ) 

11 

6 

m7 

D 

2hx2h 

6 

"8 

D 

2Hx2% 

io(r) 

All  mirrors  are  rirst  surface  aluminized  mirrors. 

(a)  Spherical  mirror. 

(b)  Diagonal  mirror (plane) . 

(c)  Distance  to  the  next  mirror. 

(d)  The  distance  from  the  exit  slit  to  is  24  inches. 

(e)  The  distances  from  M4  to  standard  mirror  and  from  standard 
mirror  to  Mg  are  both  seven  inches 

(f)  This  is  the  distance  from  Mg  to  the  detector. 


O 
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Holder  are  mounted. 

All  mirror  mounts  are  capable  of  being  adjusted  as  to 
the  height,  the  tilt,  and  the  rotation  of  the  mirror  so  the 
optics  can  be  aligned.  All  mirrors  are  cleaned  when  neces¬ 
sary  with  95%  ethanol,  distilled  water,  and  absorbent  cotton 


D.  P-E  12-C  System 

The  conventional  Perkin-Elmer  12-C  single  beam,  single 
pass,  infrared  spectrometer  with  a  chopper  in  front  of  the 
globar  source,  is  usedc  The  globar,  operated  at  about 
1100°C,  and  the  Nad  prism  give  the  desired  radiation  with 
wavelengths  ranging  from  two  to  13  microns. 

The  two  bilateral  slits  (entrance  and  exit  slits)  in 
the  P-E  Model  98  monochromator  are  each  12  millimeters  high 
which  is  reduced  to  two  millimeters  by  baffles  installed 
near  the  slits.  The  slits  are  operated  simultaneously  by 
the  slit  micrometer  which  reads  directly  in  microns  from 
zero  to  two  millimeters. 


E.  The  Detection  System 

The  detector  used  in  this  experiment  is  a  Charles  M. 
Reeder  &  Company,  Inc.,  No.  RP3-W  thermopile  with  KBr  win¬ 
dow  and  18  ohms  resistance.  The  receiver  area  is  6.0  x  0.4 
millimeters.  The  output  is  amplified  and  rectified  by  a 
Perkin-Elmer  amplifier  (Model  107)  and  monitored  on  a  Leeds 
&  Northrup  Recorder  (Model  Speedomax  W,  single  point. 


multiple  range,  multiple  speed,  strip  chart) . 

The  entire  optical  system  is  within  a  black  cloth 
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enclosure  which  prevents  dust  deposits  and  prevents  the 
detector  from  receiving  stray  light  which  tends  to  increase 
the  noise  of  the  detector  output. 

Flux  averaging  spheres,  with  two  and  three  inches  diam¬ 
eter  and  coated  inside  with  crystex  sulfur  (18),  are  fabri¬ 
cated  and  can  be  mounted  over  the  cfetector  (see  Appendix  C) . 
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VI.  CALIBRATION  OF  THE  APPARATUS 

A.  General  Features 

After  the  system  was  fabricated,  the  optics  were  proper¬ 
ly  aligned.  The  apparatus  was  calibrated  on  its  slit  width, 
wavelength,  and  linearity.  The  standard  comparison  is  then 
made  and  performance  checked. 

The  optical  alignment  of  mirrors  and  specimen  holder 
was  accomplished  by  setting  the  monochromator  so  that  the 
visible  part  of  the  spectrum  was  focused  on  the  exit  slit. 

In  place  of  the  standard  mirror  a  greased  piece  of  paper  was 
installed.  The  image  of  the  exit  slit  was  visually  focused 
at  the  center  of  the  paper  by  adjustments  of  the  mirrors, 
adjustable  screws  on  the  supporting  base  of  the  specimen 
holder  and  the  leveling  screws  of  the  specimen  holder.  Sub¬ 
sequently,  the  greased  piece  of  paper  was  installed  right 
in  front  of  the  detector.  The  image  of  the  exit  slit  was 
again  caused  to  fall  on  the  same  location  for  SPECIMEN  IN 
and  SPECIMEN  OUT  conditions.  Since  the  spherical  mirrors 
were  used  off-axis,  the  image  of  the  exit  slit  was  slightly 
astigmatic.  Optical  alignment  was  made  such  that  the  de¬ 
sired  location  was  the  circle  of  least  confusion.  The 
tangential  and  sagittal  images  were  located  very  close  to 


the  specimen  surface. 


B.  Slit  Width  Linearity 

Although  no  occasion  was  anticipated  to  compare  energies 
with  different  slit  widths,  it  was  still  necessary  to  check 
the  slit  width  linearity  to  verify  proper  operation  of  the 
monochromator  including  resolved  slit  width. 

For  optimum  alignment  the  square  root  of  the  radiative 
energy  emerging  from  the  exit  slit  should  be  proportional  to 
the  slit  width  setting  (see  Appendix  A) .  The  performance  was 
checked  at  a  wavelength  of  about  2.25/i.  At  this  wavelength 
the  globar  source  gives  the  highest  energy  and  hence,  small 
slit  widths  could  be  used.  For  each  amplifier  gain  (gains 
of-  4-0,  6-0,  9-0,  12-3),  the  slit  width  was  varied  and  the 
pen  deflections  were  recorded.  A  plot  of  the  square  root  of 
the  pen  deflection  against  the  slit  width  is  shown  in  Fig¬ 
ure  10.  The  fairly  good  linear  relation  indicates  the  good' 
alignment  of  the  monochromator  slits  and  optics. 

C.  Wavelength  Calibration 

The  drum  on  the  monochromator  has  100  divisions  for 
each  turn  with  total  of  24  turns.  The  correspondence  between 
the  drum  number  and  wavelength  is  unique  for  every  instru¬ 
ment  since  it  varies  for  each  individual  prism.  Therefore, 
it  is  necessary  to  obtain  a  wavelength  calibration  curve 
for  the  experimenter's  particular  instrument. 

The  calibration  is  in  principle  a  function  of 


Gain  12-3  Gain  9-0 


SLIT  WIDTH  (/ 1 ) 


Figure  10.  Monochromator  slit  width  Linearity 

Performance  Curve 
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temperature  of  the  prism.  The  magnitude  of  the  temperature 
effect  for  any  prism  may. readily  be  estimated  by  comparison 
of  the  variations  of  the  refractive  index  with  wavelength  and 
with  temperature.  The  temperature  immediately  adjacent  to 
the  Nacl  prism  inside  the  monochromator  has  been  checked  and 
the  temperature  was  found  to  be  controlled  at  99°F  +  1°F 
during  the  three  days'  checking.  Because  the  temperature  of 
the  prism  is  well  controlled  within  +1°F  and  the  temperature 
compensating  device  is  included  in  the  monochromator,  it  is 
expected  that  once  the  wavelength  calibration  is  established, 
there  will  be  no  change  of  this  calibration  due  to  tempera¬ 
ture  effect. 

Absorption  bands  with  known  wavelength  of  polystyrene 
(41),  H20  and  CC>2  in  atmosphere  (39)  are  used  for  the  cali¬ 
bration.  The  drum  is  driven  automatically  with  specified 
scanning  speed  to  scan  over  the  entire  spectrum.  Drum 
numbers  for  minimum  pen  deflection  are  read  precisely  by 
using  the  wavelength  marker  which  gives  a  short  line  on  the 
recorder  for  each  drum  division.  The  corresponding  wave¬ 
lengths,  drum  numbers,  together  with  the  operation  condi¬ 
tions  are  given  in  Table  2  and  plotted  in  Figure  11.  The 
corresponding  drum  numbers  for  wavelengths  from  two  microns 
to  13  microns  obtained  from  the  smoothed  calibration  curve 
in  Figure  11  is  given  in  Table  3, 


Table  2.  Wavelengths  and  Corresponding  Dram  Numbers 
Determined  From  Absorption  Bands 


(I)  Polystyrene^^ 

XjJlZ  )  Drum  Number 

2.170  1869.6 

2.577  1848.6* 

2.674  1845.2 

2.763  1841.2* 

3.269  1817.6 

3.303  1816.2* 

3.511  1806.6* 

4.225  1767.8 

4.281  1764.8* 


(II)  polystyrene^) 


Mm) 

Drum  Number 

5.138 

1712.5 

5.343 

1699.5 

5.549 

1685.0 

6.238 

1632.5 

6.692 

1596.5 

8.662 

1405.0* 

9.742 

1281.0* 

11.035 

1107.0* 

(a)  Slit  width  0.016  mm., 
rev./4  min. 

(b)  Slit  width  0.300  mm., 
rev./2  min. 

(c)  Slit  width  0.036  mm., 
rev,/4  min. 

(d)  Slit  width  0.300  mm., 
rev./2  min. 


III)  H20  in 

Atmosphere  ^C' 

tUil _ 

Drum  Number 

5.204 

1708.8* 

5.356 

1699.0 

5.466 

1690.8 

5.639 

16n8 .4* 

5.762 

1669.4 

5.822 

1665.2 

5.936 

1675.0 

5.986 

1652.8* 

6.112 

1643.2 

6.182 

1637.6 

6.339 

1625.2 

6.427 

1619.8 

6.487 

1613.2 

6.563 

1606.8* 

6.633 

1601.6 

6.709 

1595.4 

6.786 

1588.4 

6.856 

1582.2 

6.961 

1573.4 

7.044 

1565.4* 

7.165 

1554.8* 

(IV)  C02  in 

(d) 

Atmosphere' 

JlU£l _ 

Drum  Number 

13.879 

641.0* 

amplifier 

response 

2, 

drum 

speed 

amplifier 

response 

1, 

drum 

speed 

amplifier 

response 

2, 

drum 

speed 

amplifier 

response 

2, 

drum 

speed 

Data  used  to  plot  the  calibration  curve  in  Figure  11. 


Note  (1):  Settings  of  wave  drum  were  always  made  toward 
increasing  drum  number. 

Note  (2):  Polystyrene  sheet  was  mads  by  Beckman,  Inc. 
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Table  3.  Monochromator  Wavelength  Calibration. 


X(M) 

Drum  Number 

2 

1875 

3 

1830 

4 

1781 

5 

1722 

6 

1652 

7 

1565 

8 

1475 

9 

1368 

10 

1250 

11 

1120 

12 

970 

13 

800 

P .  Linearity 

Reflectivity  determination  requires  radiant  flux  meas¬ 
urements.  If  really  accurate  measurements  are  to  be  made, 
it  is  necessary  that  either  the  output  signal  from  the  in¬ 
strument  be  strictly  proportional  to  the  incident  radiant 
flux,  or  that  the  departure  from  linearity  be  known  and 
corrections  be  made.  The  linearity  of  different  stages  of 
electronics  may  be  checked  by  introducing  calibration  signals 
into  the  recorder  and  the  preamplifier  inputs  and  measuring 
the  output  signal.  The  linearity  of  the  complete  system 
(detector  and  electronics)  may  be  checked  by  introducing 
into  the  light  path  some  device  which  gives  known  trans¬ 
mittance. 
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1.  Recorder  Linearity  Calibration 
The  linearity  calibration  was  performed  by  using  a  dc 
source  (wiuh  10  ohm  output  impedance)  to  get  the  recorder 

indicator  on  100%,  90%,  - ,  and  0%  on  the  L&N  Speedomax  W 

recorder  with  rwo  millivolt  full  scale.  The  true  voltage 
value  was  read  from  the  potentiometer  and  the  procedure  was 
repeated  with  increasing  dc  source  output. 

The  calibration  curves  for  both  settings  are  shown  in 
Figure  12.  The  differences  between  the  readings  of  Recorder 
(%)  and  True  (%)  are  well  within  the  +0.3%  specified  in  the 
instrument  manual. 

Other  possible  instruments  for  readout  of  the  output 
are  digital  voltmeter  (DVM)  and  voltage-frequency  converter 
with  counter  (6).  In  the  case  of  the  counting  technique 
the  setting  precision  can  be  as  good  as  0.025%. 

Due  to  the  noise  of  the  output  signal  it  was  difficult 
to  observe  the  output  signal  on  the  displacing  screen  of  the 
above-mentioned  two  instruments  and  obtain  an  average  value. 
On  the  recorder,  however,  an  average  value  was  rather  easy 
to  obtain  from  the  output  curves  on  the  chart  paper.  Thus, 
in  this  experiment  the  recorder  was  employed. 

2.  Amplifier-Recorder  Linearity  Calibration 
The  tost  microvolts  supplied  by  the  sotting  of  the  P-E 
107  amplifier  dials  are  not  accurate  because  of  the  uncer¬ 
tainty  (~+10%)of  the  resistors  used  in  the  test  microvolt 
circuit  and  the  voltage  variation  of  the  battery.  A 
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precision  resistor  box*  and  a  2.1  volts  low  self-discharge 
rate  storage  battery  were  used  to  replace  the  dr^.  battery 
and  test  resistors  in  the  P-E  107  amplifier.  Accurately 
known  small  signal  increments  can  then  be  introduced  at  the 
input  to  the  preamplifier. 

The  Dekabox  setting  and  corresponding  True  %  are  shown 
in  Table  4.  The  calibration  results  are  shown  in  Figure  13 

Table  4.  Amplifier-Recorder  Linearity  Dekabox  Setting 


R(Kfi) 


True  % 


10 

12.5 

16 

20 

25 

40 

50 

100 

open 


100 

80 

62.5 

50 

40 

25 

20 

10 

0 


for  both  increasing  and  decreasing  voltage  settings.  The 
apparent  (%) -True  (%)  are  within  j4D.2%. 

3.  Detector-Amplifier-Recorder  Linearity  Calibration 
The  final  calibration  required  is  that  of  the  linoaritv 
of  the  complete  system  (detector  and  electronics).  This 
calibration  can  be  made  for  the  thermopile  detector  by  using 


★ 

Model  DD  655  Dekabox  decade  resistors  are  double  coaxial- 
dial  precision  resistance  elements  manufactured  by  Electro 
Scientific  Industries,  Inc.,  Portland,  Oregon.  The  accuracy 
for  lO^  ohms  to  10^  ohms  is  0.01%. 
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a  set  of  sector-disc  attenuators  similar  to  the  set  described 
by  Harrison,  et  al.,  (28).  The  discs  are  the  same  as  used 
by  Zipin  (J")  with  the  detailed  dimensions  and  true  trans¬ 
missions  given  in  Table  5. 

Table  5.  True  Transmission  of  the  Sector-Disc  Attenuatbrs 


Nominal 
Trans¬ 
mission 
of  Disc(%) 

Diameter 
(in. ) 

Sum  of 

Notches 

Angle 

l  a 

True 
Trans¬ 
mission 
T(%) _ 

No.  of 
Notches 

Uncertainty 
A  (T(%)] 

5 

10 

22°27 1 

6.24 

4 

+0.0185 

12.5 

9 

43°46 ' 

12.16 

8 

+0.0370 

25 

9 

90°9  1 

25.58 

8 

+0.0370 

50 

9 

185°42 1 

51.58 

8 

+0.0370 

75 

9 

271°10 1 

75.32 

8 

+0.0370 

To  perform  the  calibration,  the  specimen  holder  is  set 
in  the  SPECIMEN  OUT  position  and  the  amplifier  gain  is  ad¬ 
justed  to  bring  the  recorder  deflection  to  near  full  scale. 
The  discs  are  put  in  the  optical  path,  between  M3  and  M4,  in 
turn  and  driven  with  a  motor*  at  a  high  enough  speed 
(~1400  rpm)  so  that  no  coupling  is  observed  between  the 
sectored  disc  and  the  chopper.  The  discs  are  not  put  right 
in  front  cf  detector  so  there  will  be  no  thermal  radiation 
falling  on  the  detector  directly  due  to  the  presence  of  the 
motor  and  disc  blades  to  increase  the  noise.  The  recorder 
deflection  due  to  each  of  the  discs  rotating  in  the  optical 


*Manufactured  by  precision  Scientific  Co.,  Chicago,  Illinois. 
The  motor  is  turning  counterclockwise  and  the  chopper  blade 
is  turning  clockwise  when  observed  along  the  optical  path. 
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path  is  compared  with  the  100%  reading  to  calculate  the 
apparent  transmission.  They  are  plotted  as  [t  apparent  (%)  - 
T  true  (%)]  vs  (T  apparent)  in  Figure  14.  The  non-linearity 
is  about  +0.6%. 

The  operation  principle  of  the  attenuator  is  that  the 
radiation  energy  arriving  at  the  thermopile  detector  is 
blocked  by  the  rotating  blades  of  the  attenuator  part  of  the 
time.  The  portion  of  tne  unblocking  time  determines  the 
transmission.  There  is  •inherent  uncertainty  involved  in 
using  the  attenuator  as  a  device  to  check  the  linearity  be¬ 
cause  of  the  slower  response  time  (comparative  co  electronics) 
of  the  thermopile.  That  is,  the  real  transmission  caused  by 
the  attenuator  may  never  be  the  same  as  the  true  transmission 
value  given  in  Table  5,  no  matter  how  small  the  uncertainty 
of  the  angles  of  the  blades  can  be  determined.  It  may  seem 
possible  at  the  first  thought  to  make  these  two  transmissions 
identical  by  using  a  slower  speed  of  the  attentuator  motor. 

But  in  this  way,  the  coupling  effect  with  the  chopping  blades 
in  front  of  the  source  increases  the  noise  of  the  output  sig¬ 
nal  and  thus  increases  the  uncertainty  of  the  linearity 
determination. 

The  other  devices  not  involving  rotating  discs,  such  as 
standard  filters,  although  not  handicapped  by  the  uncertainty 
mentioned  above,  still  yield  systematic  errors  (3).  still 
another  device,  the  line  screen  (or  mesh  screen)  does  not 
give  rise  to  most  of  the  systematic  errors,  but  ir  is 


(%) 3nyi 


(%)  iusjoddv 


Figure  14 *  Linearity  Characteristic  of  the  Detector  System. 
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very  sensitive  to  the  convergence  angle  of  the  optical  beam. 
Thus,  by  these  three  above-mentioned  devices,  it  is  diffi¬ 
cult  to  determine  linearity  to  better  than  1%  or  2%.  The 
most  promising  device  is  a  system  employing  three  optically 
excellent,  high-extinction  polarizers* (3) .  The  common 
sources  of  the  systematic  error  are  eliminated  and  the 
linearity  can  be  determined  to  better  than  0.1%. 


Stray  energy  may  often  cause  important  errors  in  spec¬ 
trometry.  The  stray  energy  emerging  from  the  exit  slit  along 
with  the  desired  spectral  energy  may  be  a  general  admixture 
of  energy  of  all  wavelengths  emitted  by  the  source  and  trans¬ 
mitted  by  the  spectrometer. 

Stray  energy  was  checked  by  using  numbers  210  and  240 
infrared  filters  purchased  from  Eastman  Kodak  Company.  Both 
were  two  inches  in  diameter  x  0.04  inches  thick.  No.  210 
and  no.  240  infrared  filters  are  filters  without  polystyrene 
protective  coating  with  1M  and  4M  cut-on  wavelength**, 
respectively.  The  infrared  filter  was  placed  in  front  of  the 
entrance  slit  of  the  monochromator  during  the  measurements,. 
The  No.  210  filter  was  used  for  the  wavelength  range  two  to 
five  microns;  while  the  No,  240  filter  was  used  in  the  seven 


to  13  micron  range. 


*  An  instrument  employing  such  optical  elemer.cs  is  called 
the  Linearity/Transmittance  Standards  (L/TS)  now  in  pro¬ 
duction  by  Technometrics,  Inc.,  West  Lafayette, Indiana. 

**The  cut-on  wavelength  is  defined  as  the  wavelength  at  which 
the  transmission  reaches  3%. 


Reflectivity  of  aluminum  specimen  No.  A^-Bl  was  meas¬ 
ured  with  and  without  filter  from  two  to  13  microns.  No 
trend  of  difference  can  be  detected  from  the  results  (see 
Figure  15)  measured  with  these  two  conditions.  Since 
rhodium  has  lower  reflectivity  value,  the  reflectivity  of 
rhodium  specimen  No.  1  was  thus  measured  with  and  without 
filter  in  the  wavelength  range  seven  to  13  microns  (see 
Figure  16)  in  order  to  check  whether  there  exists  any  differ 
ence 

The  results  of  the  reflectivity  measurements  with  and 
without  the  filter  all  fall  within  the  experimental  uncer¬ 
tainty,  indicating  that  the  stray  energy  emerging  from  the 
exit  slit  is  negligible.  Similar  experience  is  reported  by 
Edwards  and  devolo  (19)  with  this  type  of  apparatus. 


Figure  15.  Stray  Energy  Determination 

(A-l*  Reflectivity) 
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VII.  EXPERIMENTAL  PROCEDURE 


Upon  completion  of  the  alignment  and  calibration  of  the 
apparatus,  the  standards  comparison  and  performance  check 
were  then  made.  After  the  surface  preparation  of  the  speci¬ 
mens,  the  measurements  on  the  specimens  were  performed.  The 
standards  comparison  and  surface  preparation  will  be  dis¬ 
cussed  in  Chapters  VIII  and  IX. 

The  whole  electronic  instrumentation  system  (except  the 
globar  source)  was  kept  on  during  the  entire  measuring  period 
to  avoid  any  possibility  of  transient  effects.  After  turning 
on  the  globar  source,  there  was  usually  a  wait  of  approxi¬ 
mately  one  hour  before  it  reached  steady  state.  The  stability 
of  the  globar  source  was  checked  by  observing  minimum  pen  de¬ 
flection  on  the  recorder  for  a  ten  minute  period. 

First,  the  specimen  and  two  small  standard  mirrors  were 
mounted  on  the  specimen  holder  for  the  orientation  AB-Rl . 

The  "zero"  position  on  the  recorder,  i.e.  with  shutter  in 
front  of  the  entrance  slit  closed  and  detector  uncovered, 
was  then  set  by  adjusting  the  BALANCE  dial  on  the  amplifier. 

A  two  millivolt  span  was  used  on  the  recorder  and  one  milli¬ 
volt  was  suppressed. 

The  desired  wavelength  drum  setting,  slit  width,  and 
amplifier  gain  were  set,  and  the  pen  deflections  for  SPECIMEN 
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IN  and  SPECIMEN  OUT  positions,  as  described  in  Chapter  V, 
were  recorded  for  each  wavelength.  Measurements  were  per¬ 
formed  for  wavelengths  2,3,4,5,7,9,11,  and  13  microns.  The 
slit  width  setting  varied  from  0.080  millimeters  at  two 
micron  wavelength  to  1.58  millimeters  at  13  microns  corre¬ 
sponding  to  0.035  micron  and  0.28  micron  resolved  bandwidth, 
respectively  (see  Appendix  B) .  After  all  three  specimens 
were  measured  for  one  block  orientation,  the  same  procedure 
was  carried  out  for  the  other  three  orientations.  Every 
specimen  was  measured  within  eight  hours  after  surface  pre¬ 
paration. 

The  flux  averaging  spheres,  two  and  three  inch  diameters 
with  crystex  sulfur  coating,  were  tried  over  the  detector. 

The  output  was  too  small  especially  at  the  longer  wavelengths 
(see  Appendix  C) .  Therefore,  it  is  hard  to  obtain  an  ade¬ 
quate  value  to  calculate  the  radiant  flux  ratio  precisely. 
Thus,  all  the  measurements  were  made  without  the  flux  averag¬ 
ing  sphere.  A  source  of  error  arising  from  spatial  sensitiv¬ 
ity  of  the  detector  cannot  be  avoided.  Hence,  the  determin¬ 
ation  of  measurement  accuracy  and  precision  is  strongly  re¬ 
quired  and  will  be  discussed  in  the  following  section  "Stan¬ 
dards  Comparison  and  Performance  Check." 

The  data  reduction  procedure  is  as  follows.  For  each 
specimen  the  reflectivity  was  measured  as  a  function  of  wave¬ 
length.  For  every  orientation  and  at  each  wavelength  the  re¬ 
flected  energy  of  SPECIMEN  IN  and  SPECIMEN  OUT  were  recorded 
to  calculate  the  normal  spectral  reflectivity.  Since  the 
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same  amplifier  gain  was  used  for  SPECIMEN  IN  and  SPECIMEN 
OUT  positions,  no  determination  of  instrument  response  con¬ 
stant  was  needed. 

In  this  experiment  the  recorder  was  adjusted  to  zero 
when  the  shutter  in  front  of  the  entrance  slit  was  closed. 
Thua ,  the  dc  signal  due  to  the  unwanted,  scattered  radiation 
reaching  the  detector  and  stray  electrical  signals  were  al¬ 
ready  excluded  from  the  recorder  readings.  For  every  orien¬ 
tation,  the  normal  spectral  reflectivity,  R^,  is  calculated 
by 

Ri  “  <VD2)!S 

where  and  d2  are  recorder  pen  deflections  corresponding 
to  SPECIMEN  IN  and  SPECIMEN  OUT  positions. 

The  final  value  for  normal  spectral  reflectivity  is  cal¬ 
culated  from  data  obtained  from  four  different  orientations 
R  =  (R1R2R3R4)4 


VIII.  STANDARDS  COMPARISON  AND  PERFORMANCE  CHECK 


Reproducibility  of  the  measurements  was  checked  by  com¬ 
paring  results  of  the  first  run  and  second  run  of  the  reflec¬ 
tivity  measurements  on  opaque,  uncharacterized  rhodium  speci¬ 
men  No.  1*  Rhodium  was  selected  because  its  reflectivity 
characteristics  are  very  stable  in  air  and  at  room  tempera¬ 
ture.  The  specimen  was  electroplated  rhodium  on  flat  glass 
prepared  by  Evaporated  Metal  Films  Corporation,  Ithaca,  New 
York.  The  results  of  two  runs  are  well  within  0.002  reflec¬ 
tivity  unit  of  each  other  (see  Figure  17).  Results  are  also 
comparable  to  the  data  reported  by  Hass  (29). 

The  accuracy  of  the  measurement  technique  was  evaluated 
as  follows..  Three  quartz  flats  were  coated  with  high-vacuum 
evaporated  aluminum*  under  well  controlled  conditions  which 
yielded  specimens  of  reference  grade.  The  reflectivity  of 
these  three  specimens  v/ere  measured  in  this  study  two  weeks 
and  four  weeks  after  evaporation.  The  average  values  are  shown 
in  Figure  18.  A  fourth  specimen,  prepared  at  the  same  time 
as  the  test  specimens  referred  to  above  was  measured  immedi¬ 
ately  after  evaporation  at  the  Michelson  Laboratory,  and 

Prepared  by  Dr.  H.  E.  Bennett,  at  the  Michelson  Laboratory, 
Naval  weapons  Center,  China  Lake,  California. 


1st  run 


Figure  17.  Reproducibility  Evaluation  From  Measurements  on 

The  Reflectivity  of  Rhodium. 
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Figure  18.  Accuracy  Evaluation  from  Measurements  on 

the  Reflectivity  of  Aluminum. 
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data  obtained  are  plotted  on  Figure  18,  together  with  earlier 
literature  values  (5)  on  fresh  and  aged  high  vacuum  evapor¬ 
ated,  aluminized  specimens.  Generally,  the  reflectivity 
measurements  for  the  condition  "four  weeks  after  evaporation" 
are  within  0.001  reflectivity  unit  from  those  of  Bennett 
(5)  stated  to  have  an  accuracy  of  0.001  unit.  However,  the 
value  of  reflectivity  obtained  in  this  experiment  is  system¬ 
atically  higher  at  the  longer  wavelengths.  This  cannot  be 
attributed  to  the  stray  energy  emerging  from  the  exit  slit 
since  it  is  negligible  as  demonstrated  in  Chapter  VI,  "Cali¬ 
bration  of  the  Apparatus." 


f 
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IX.  SPECIMEN  PREPARATION  AND  SURFACE  CHARACTERIZATION 


A.  Specimen  Materials 

Copper,  zinc,  and  nickel  rods  1.25  inches  in  diameter 
and  two  inches  long  were  purchased  from  Material  Research 
Corporation,  Orangeburg,  New  York.  The  purities,  as  reported 
by  the  manufacturer,  are  given  in  Table  6. 


Table  6.  Purity  of  the  Specimens 


Material3 

Purity (%) 

Grade 

Total  Impurity (ppm) c 

Copper 

99.999 

MARZb 

<33 

Zinc 

99.999 

MARZb 

<  6 

Nickel 

99.99 

VP 

<123 

a.  Vacuum  melted  and  cold  swaged  to  size,  certificated. 

b.  zone  refined. 

c.  According  to  supplier's  typical  emission  spectrographic 
analysis. 

B.  Surface  Preparation 

All  rods  were  first  longitudinally  milled  flat  on  two 
sides  with  a  resulting  face  width  of  1.125  inches  in  order  to 
fit  into  the  specimen  holder.  Four  specimens  of  0,25  inch 
thickness  each  were  cut  from  each  rod.  For  each  metal,  one 
specimen  was  first  polished  by  using  metallographic  methods 
to  test  the  effectiveness  of  the  polishing  procedure.  If 
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the  specimen  surface  thus  prepared  was  smooth  and  bright,  the 
remaining  three  specimens  were  prepared  by  using  the  same 
procedure.  The  procedure  used  was  as  follows. 

Specimens  of  0.25  inch  thickness  were  cut  off  from  the 
rods  with  a  Felker  (Torrance,  California)  Di-Met  Model  HR 
cut-off  machine.  A  resinoid  bonded,  metal  oxide  cut-off 
wheel,  A70-R8-B1,  8  inch  diameter  and  0.0625  inch  thick,  was 
used  for  the  nickel  rod  and  a  Electroion  silico  carbide  cut¬ 
off  wheel,  C80-010-B1,  8  inch  diameter  and  0.1875  inch  thick, 
was  used  for  both  the  copper  and  zinc  rods.  Different  cut¬ 
ting  techniques  were  required  as  nickel  has  a  BrinelD  hard¬ 
ness,  soft  condition  (500  kg),  of  70  while  that  for  copper 
and  zinc  are  30  and  35,  respectively.  The  cut-off  wheels  are 
manufactured  by  Simonds  Abrasive  Company,  Philadelphia, 
Pennsylvania  and  were  operated  at  3600  rpm.  Sufficient  cool¬ 
ing  water  was  used  during  the  cutting  to  prevent  any  heat 
tempering  to  avoid  alteration  of  the  microstructure  and  hard¬ 
ness  of  the  specimen.  During  cutting,  a  reasonable  amount  of 
wheel  pressure  was  used  so  the  abrasive  grains  would  "bite" 
and  remove  the  metal  in  fine  chips,  and  at  the  same  time,  • 
cause  the  wheel  bond  to  break  down  and  present  new  abrrsive 
particles  to  .the  cutting  surface. 

The  metallographic  polishing  consists  of  coarse  and  fine 
grindings,  rough  polishing,  and  final  polishing.  Nos.  240, 
320,  400,  and  600  grit  silicon  carbide  papers  were  used,  and 
water  was  provided  during  the  coarse  and  fine  grinding 
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processes.  Each  step  was  concluded  when  the  surface 
appeared  uniformly  dull  over  its  entire  area  using  visual 
examination. 

Rough  polishirg  was  performed  by  placing  the  specimen 
on  a  nylon  cloth  over  the  polishing  wheel  with  6M  Buehler 
Metadi  diamond  as  the  polishing  compound  and  Buehler  No. 
60-3250  AB  lapping  oil  as  the  cooling  liquid.  The  aim  of 
this  stage  was  to  eliminate  any  disturbed  surface  layers  and 
to  prepare  an  overall  surface  flatness.  After  polishing, 
the  specimen  was  washed  first  with  trichloroethylene,  then 
with  ethanol,  and  finally  with  tap  water. 

Gamma  form  of  aluminum  oxide  with  distilled  water  was 
used  on  the  Buehler  microcioth  for  the  final  polishing  stage 
Microcloth  is  a  napped  cloth  which  is  compressible  and  tends 
to  conform  to  the  surface  even  under  the  slightest  pressure. 
The  cloth  action  coupled  with  the  abrasive  action  has  a 
greater  removal  rate  on  the  softer  phases  than  on  the  harder 
phases. 


C.  Surface  Examination 

The  requirements  of  the  surface  preparation  are  that  the 
resulting  surface  should  be  flat  to  minimize  the  error  aris¬ 
ing  from  spatial  sensitivity  of  the  detector  and  smooth  to 
reduce  the  diffuse  component  which  gives  rise  to  experimen¬ 
tal  error. 

Great  care  has  been  taken  during  the  grinding  and  polish' 


ing  processes  in  order  to  prepare  a  flat  surface  on  the 


specimen.  Those  specimens  with  good  edge  preservation  were 
assumed  to  be  sufficiently  flat  and  were  selected  for  the 
reflectivity  measurements. 

If  the  surface  is  rough,  the  reflected  flux  from  the 
specimen  which  is  within  the  viewing  solid  angle  of  the 
detection  system  is  reduced  by  virtue  of  a  larger  hemispher¬ 
ical  diffuse  component.  The  effect  of  roughness  is  to  re¬ 
distribute  the  reflected  flux  so  that  the  rough  surface  has 
an  increased  diffuse  component  in  addition  to  the  predomin¬ 
ant  specular  component.  For  the  roughness  range  being  con¬ 
sidered  here,  the  effect  is  one  of  redistributing  the  flux- 
and  not  increasing  the  absorptivity  (2). 

The  normal  spectral  reflectivity  is  a  very  sensitive 
indicator  of  roughness  when  the  wavelength  is  much  larger 
compared  with  the  roughness.  The  relation  is  given  as 

R/Ro  =  exp[- (4ffa ) 2 /A2] 

where  R  is  the  specular  reflectance  of  the  rough  surface 
with  surface  roughness  o  and  Rq  that  of  a  perfectly  smooth 
surface  of  the  same  material.  From  this  relation  the  error 
caused  by  surface  roughness  in  reflectance  measurement  was 
estimated  (2).  For  the  wavelength  range  of  this  experiment, 
the  surface  roughness  should  '  e  less  than  about  0.02  micron 
in  order  to  give  less  than  0.1?4  error  in  reflectance.  The 
final  stage  of  polishing  of  this  experiment  was  performed 
by  using  gamma  alumina  slurry  which  has  particle  size  of 
0.05  micron.  Thus,  the  requirement  can  be  reasonably 
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achieved. 

In  this  range  of  roughness,  the  conventional  technique 
such  as  mechanical  profilometer  is  inaccurate  and  it  is  also 
destructive  to  the  surface  measured.  The  surface  roughness 
can  be  examined  by  using  the  Thermal  comparator*.  The  Ther¬ 
mal  Comparator,  a  tool  for  rapid  measurement  of  thermal  con¬ 
ductivity,  is  also  quite  useful  for  examination  of  surface 
characteristics  and  for  qualitative  determinations  of  film 
thickness,  roughness,  oxidation,  etc. 

The  Thermal  Comparator  consists  of  a  thermocouple 
mounted  very  close  to  the  surface  of  a  fine  pointed  probe 
which  in  turn  is  mounted  in  a  large  heat  source  maintained 
at  some  excess  temperature  (15-30°C)  above  the  ambient.  When 
the  probe  is  touched  to  a  specimen  surface,  the  tip  at  which 
the  thermocouple  (differentially  connected  with  the  main 
source)  is  located,  quickly  reaches  an  intermediate  tempera¬ 
ture  between  ambient  and  that  of  the  source.  If  the  speci¬ 
men  is  a  smooth,  homogeneous  solid,  the  ratio  of  the  thermo¬ 
couple  EMF  from  the  specimen  when  compared  to  a  similarly  pre 
pared  surface  of  a  standard  reference  material  will  yield  the 
thermal  conductivity.  If  the  surfaces  are  roughened,  the 
point  probe  will  contact  the  asperities  rather  than  the  full 
surface  resulting  in  a  lower  EMF  response  and  hence  giving  a 
qualitative  measure  of  roughness.  The  Thermal  Comparator  can 

•fa 

Manufactured  by  the  Technometrics,  Inc.,  West  Lafayette, 

Indiana. 
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also  be  used  to  rapidly  determine  the  similarity  of  surface 
preparation  for  a  collection  of  specimens  before  reflectance 
measurements . 

To  determine  the  influence  of  the  polishing  technique 
on  roughness,  measurements  of  the  Thermal  Comparator  response 
were  made  following  the  various  stages  of  polishing.  Polish¬ 
ing  continued  until  all  specimens  indicated  the  same  EMF 
(qualitatively  having  the  same  minimum  roughness). 

Zinc  specimens  were  purposely  prepared  with  different 
resulting  surface  roughness.  The  results  of  the  reflectivity 
measurement  and  the  Thermal  Comparator  EMF  response  of  each 
specimen  are  shown  in  Figure  19.  As  reasoned  above,  the 
zinc  specimen  (No. 2)  with  highest  EMF  response  gives  the 
highest  reflectivity. 

From  the  reflectivity  data  given  in  Figure  19,  the  sur¬ 
face  roughness  of  specimens  Nos.  1  and  3  (assuming  specimen 
No.  2  as  a  reference)  by  using  the  relation  of  the  reflec¬ 
tance  and  roughness,  at  four  and  eight  microns,  gave  the 
roughness  value  as  0.02-0.04  micron  for  specimen  No.  1  and 
0.06-0.10  micron  for  specimen  No.  3. 


8 


X0  EXPERIMENTAL  RESULTS  AND  THEIR  COMPARISON 

WITH  THEORY 


Measurements  of  normal  spectral  reflectivity  were  made 
on  metailographically  polished  copper,  zinc,  and  nickel,  and 
aged,  high  vacuum  evaporated  aluminum  specimens.  Three  speci¬ 
mens  for  each  polished  material,  which  visual  examination 
showed  good  edge  preservation,  were  selected  for  the  measure¬ 
ment.  Figures  20,  23,  and  26  show  the  measurement  results  of 
copper,  zinc,  and  nickel.  The  comparison  between  the  present 
data  and  literature  values  are  presented  in  Figures  21,24,  and 
27.  Data  on  evaporated  aluminum  and  their  comparison  with 
that  of  Bennett'  have  been  shown  in  Figure  18  in  the  previ¬ 
ous  chapter.  The  comparisons  of  measurements  and  the  predic¬ 
tions  of  theoretical  models  are  given  in  Figures  22,25,28, 
and  29.  They  are  plotted  as  absorptivity,  i.e.,  1-R,  vs 
wavelength  in  full  logarithmic  coordinates.  The  physical 
parameters  used  for  the  theoretical  calculations  are  listed 
in  Table  7. 


Copper  -  Average  reflectivity  values  of  specimen  Nos.  1, 
2,  and  4  are  plotted  in  Figure  20,  Generally,  deviations  of 
the  experimental  data  from  a  smooth  curve  are  rather  small. 
The  scattering  of  data  at  three  and  five  microns  is  due  to 
the  atmosphere  absorption  giving  rise  to  a  higher  noise  level 
of  the  output  signal  and  higher  uncertainty.  The  present 
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Figure  21.  Normal  Spectral  Reflectivity  of  Copper  and 
Its  Comparison  with  Published  values. 
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results  are  comparable  with  Bennett's  values  (Figure  21). 

Theoretical  predictions  and  their  comparison  with  ex¬ 
perimental  data  are  presented  in  Figure  22.  Prediction  of 
absorptivity  using  the  theory  of  anomalous  skin  effect  has 
been  made.  When  diffuse  electronic  reflection,  one  conduc¬ 
tion  electron  per  atom,  and  bulk  dc  conductivity  axe  assumed, 
the  prediction  is  lower  than  the  present  data.  No  further 
attempt  was  made  to  fit  the  experimental  data  with  the 
theoretical  equation  derived  from  the  theory  of  anomalous 
skin  effect.  The  predicted  curve  by  using  quantum  theory  is 
within  the  uncertainty  of  the  experimental  values.  Absorp¬ 
tivity  calculated  by  using  Hagen-Rubens  relation  is  presented 
to  show  the  trend  expected  toward  long  wavelengths.  Assuming 
1.35  conduction  electrons  per  atom  (from  Hall  effect  measure¬ 
ment)  (34),  the  simple  Drude  theory  (4)  yields  that,with  0.5 
bulk  dc  electrical  conductivity,  best  agreement  can  be  ob¬ 
tained  between  the  experimental  data  and  the  theoretical  pre¬ 
diction  at  13  microns.  These  two  parameters  are  used  to  pre¬ 
dict  the  absorptivity  from  5  microns  to  longer  wavelengths. 

The  reason  for  the  discrepancy  between  theory  and  experimen¬ 
tal  values  at  the  shorter  wavelength  end  is  because  in  the 
neighborhood  of  the  absorption  edge  and  at  higher  frequencies, 

the  damping  coefficient,  which  is  assumed  to  be  constant, 

/ 

becomes  frequency  dependent  (4). 

Zinc  -  Data  of  specimen  No.  2  are  presented  in  Figure 
23.  Specimen  No.  2  has  the  most  smooth  surface  among  three 
specimens  prepared.  The  experimental  curve  is  very  smooth 
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Figure  23.  Normal  Spectral  Reflectivity  of  zinc 


NORMAL  SPECTRAL  REFLECTIVITY 


85 


O 

□  □ 


Reflectivity  For  Mechanically 
Polished  zinc 

•  This  Experiment  (7°, 7°) 

O  Edwards, et  al(19)  (2tt,15°) 

A  Coblent z (10)  (15°, 15°) 

□  Dunkle,et  al(16)  (2ff,5°) 

O  Coblentz (9)  (25°, 25°) 

V  Hodgson (30) 


0.92 


0.90 


0.8  8  r 


0.86 


0.8  4  • - *■ 

2  4 


i  i  i 

6  8  10 

WAVELENGTH  U) 


12 


Figure  24.  Normal  Spectral  Reflectivity  of  zinc  and 
Its  Comparison  With  Published  Values. 
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Figure  25.  Absorptivity  of  zinc  and  Its  comparison  with  Theory. 
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except  the  lower  value  at  11  microns  which  is  probably  due  to 
the  experimental  error.  When  compared  with  literature  values 
(Figure  24),  the  present  results  are  among  the  highest. 

Theoretical  predictions  and  their  comparison  with  ex¬ 
perimental  data  are  presented  in  Figure  25.  The  predictions 
made  bv  the  theory  of  anomalous  skin  effect  and  quantum 
theory  and  their  comparisons  with  experimental  data  are 
similar  to  that  of  the  copper  specimens.  The  experimental 
value  has  the  similar  slope  with  the  curve  calculated  from 
the  Hagen-Rubens  relation.  Assuming  2.8  conduction  electrons 
per  atom  (from  Hall  Effect  measurement) (34) ,  the  simple  Drude 
theory  (4)  yields  that  with  bulk  dc  electrical  conductivity, 
best  agreement  can  be  obtained  between  the  experimental  data 
and  the  theoretical  prediction  at  13  microns.  These  two 
parameters  are  used  to  predict  the  absorptivity  from  five 
microns  to  longer  wavelengths. 

Nickel  -  The  reflectivity  data  for  the  three  nickel 
specimens  are  shown  in  Figure  26.  Because  nickel  has  a 
rather  high  hardness,  specimens  are  easily  polished  and  yield 
a  smooth  and  bright  surface.  Nickel  surface  is  not  oxidized 
in  air  at  room  temperature.  These  may  explain  why  the  data 
of  all  three  specimens  agree  so  closely  with  each  other. 

Theoretical  predictions  and  their  comparison  with  experi 
mental  data  are  presented  in  Figures  27  and  28.  Prediction 
of  absorptivity  using  the  theory  of  anomalous  skin  effect  has 
been  made.  When  diffuse  electronic  reflection, one  conduction 
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Figure  26.  Normal  Spectral  Reflectivity  of  Nickel 
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Figure  27.  Normal  Spectral  Reflectivity  of  Nicfcel  and 
Its  comparison  with  Published  Values. 
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electron  per  atoir.,  and  bulk  dc  conductivity  are  assumed,  the 
prediction  is  lower  than  the  present  data.  No  calculate  is 
hiade  based  on  quantum  theory.  Nickel  is  a  ferromagnet: r  ma¬ 
terial  and  it  is  rather  difficult  to  estimate  its  physical 
parameters  from  the  two  simple  optical  constants  and  wave¬ 
length  relations.  Assuming  0o6  conduction  electron  per 
atom  (37),  the  simple  Drude  theory  (4)  yields  that  with  0.6 
bulk  dc  electrical  conductivity,  best  agreement  can  be  ob¬ 
tained  between  the  experimental  data  and  the  theoretical  pre¬ 
diction  at  13  microns.  These  two  parameters  are  used  to 
predict  the  absorptivity  from  about  10  microns  to  longer 
wavelengths. 

Aluminum  -  Experimental  data  on  aged,  high  vacuum 
evaporated  aluminum  (Figure  18)  are  comparable  with  litera¬ 
ture  values  (5)  of  specimen  prepared  under  the  same  condition. 

Theoretical  predictions  and  their  comparison  with  experi¬ 
mental  data  are  presented  in  Figure  29.  The  predictions  made 
by  the  theory  of  anomalous  skin  effect  and  quantum  theory  and 
their  comparisons  with  experimental  data  are  similar  to  that 
of  the  copper  specimens.  Absorptivity  calculated  by  using 
Hagen-Rubens  relation  is  presented  to  show  the  trend  expected 
toward  long  wavelengths.  Bennett  et  al  (5)  have  shown  the 
validity  of  the  Drude  theory  for  ultra-high  vacuum  evaporated 
aluminum  film  in  the  infrared  range.  Their  calculation  was 
based  on  bulk  dc  conductivity  and  2.6  conduction  electrons  per 
atom  of  aluminum.  Assuming  2.6  conduction  electrons  per 
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Figure  29.  Absorptivity  of  Aluminum  and  Its  comparison  with  Theory. 
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atom,  the  simple  Drude  theory  (4)  yields  that  with  0.5  dc 
electrical  conductivity,  best  agreement  can  be  achieved  be¬ 
tween  the  present  experimental  data  and  the  theoretical  pre¬ 
diction  at  13  microns.  These  two  parameters  are  used  to  pre¬ 
dict  the  absorptivity  in  the  infrared  range. 

Smooth  curves  based  upon  present  measurement  and  theory 
are  drawn  for  copper,  zinc,  nickel,  and  aluminum.  They  are 
shown  in  Figures  30  and  31.  Absorptivity  data  extracted 
from  these  smooth  curves  are  presented  in  Table  8,  and  are 
considered  as  the  most  probable  absorptivity  values. 


Table  8.  Most  Probable  Absorptivity  Values  Based 

Upon  Measurement  and  Theory* 


Mm) 

Cu 

zn 

Ni 

At 

2 

0.0190 

0.0560 

0.1550 

0.0305 

4 

0„ 0140 

0.0342 

0.0600 

0.0236 

6 

0.0126 

0.0270 

0.0385 

0.0210 

8 

0.0118 

0.0235 

0.0302 

0.0193 

10 

0.0114 

0.0215 

0.0260 

0.0182 

15 

0.0107 

0.0189 

0.0220 

0.0161 

20 

0.0101 

0.0175 

0.0205 

0.0148 

25 

0.0096 

0.0166 

0.0192 

0.0138 

30 

0.0092 

0.0156 

0.0180 

0.0130 

35 

0.0087 

0.0147 

0.0172 

0.0122 

40 

0.0083 

0.0140 

0.0164 

0.0115 

50 

0.0073 

0.0128 

0.0150 

0.0105 

60 

0.0068 

C . 0118 

0.0140 

0.0096 

70 

0.0062 

0.0110 

0.0132 

0.0090 

80 

0.0058 

0.0103 

0.0124 

0.0084 

90 

0.0054 

0.0098 

0.0117 

0.0079 

100 

0.0051 

0.0093 

0.0111 

0.0074 

*  .  . 
Preliminary 

estimate; 

uncertainty 

+0.002  unit. 

Figure  30.  Most  Probable  Absorptivity  Curves  of  Cu  and  A-^  Based 

Upon  Measurement  and  Theory. 
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XI.  SUMMARY  AND  CONCLUSIONS 

An  apparatus  was  fabricated  to  measure  normal  spectral 
reflectivity.  The  multiple  reflection  technique  utilized 
in  this  investigation  allows  the  reflectivity  measurement 
of  infrared  radir.tion  to  have  the  accuracy  of  about  +0.002 
reflectivity  unit. 

Measurements  of  normal  spectral  reflectivity  were  made 
on  the  pure  metallic  elements:  metallographically  polished 
copper,  zinc,  and  nickel  and  aged,  high  vacuum  deposited 
aluminum.  Special  attention  was  given  to  the  performance 
evaluation  of  the  reflectometer  and  bulk  specimen  prepara¬ 
tion. 

The  theoretical  portion  of  this  study  has  included  the 
review  of  the  theoretical  models  in  the  free-electron  fre¬ 
quency  range  and  their  applicability.  It  has  been  found 
that  absorptivity  predicted  by  the  theory  of  anomalous  skin 
effect,  assuming  diffuse  electronic  reflection,  one  conduc¬ 
tion  electron  per  atom  and  bulk  conductivity,  is  lower  than 
that  of  the  present  experiment  <,  No  attempts  were  made  to 
fit  the  data  using  this  theory.  On  the  other  hand,  it  is 
seen  that  the  prediction  made  by  simple  quantum  theory  fits 
well  with  the  experimental  results.  Simple  Drude  Theory  and 
Hagen-Rubens  relation  were  used  to  predict  absorptivity 
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values  in  the  longer  wavelength  range  not  measured  in  this 
experiment.  Most  probable  absorptivity  values  of  copper, 
zinc,  nickel,  and  aluminum  from  two  to  100  microns  were 
obtained,  based  on  measurement  and  theory. 
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APPENDIX  A 


MONOCHROMATOR  SLIT  WIDTH  LINEARITY 


In  a  single  pass  monochromator,  the  radiation  beam, 
after  passing  through  the  entrance  slit,  goes  through  the 
"collimator"  to  the  prism  and  then  from  the  prism  to  the 
"telescope"  and  to  the  exit  slit,  (see  Figure  A-l). 

The  spectral  range  of  radiation  falling  on  the  exit 
slit  of  width  A  X2  in  the  focal  plane  of  the  "telescope"  is; 


A  \  AX 
A  =  TT  * 


(A-l) 


The  radiant  flux,  B,  accepted  by  an  entrance  slit  with 
A  X-^  and  height  t  is  such  that 


3  & 


<•)[(**•*  *) 


(A-2) 


Where  E-^’AXis  the  power  emitted  by  the  source  per  unit  area 
per  unit  solid  angle  in  the  wavelength  range  X  to  X+A  X, 
Aberrations  of  mirrors  and  curvature  of  the  spectral  images 
are  second  order  effects  and  may  be  safely  ignored;  losses 
by  reflection  polarization  and  absorption  are  neglected  here. 

Inserting  (A-l)  into  (A-2)  and  setting  f^-f2=f  and 
A  X-^  =  A  X2  =  s  conditions  for  the  P-E  Model  98  monochromator, 
yield 
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O 


3~[^Ex(ax/de)](s>2  <A-3> 

The  parameters  l,  D,  and  f  era  fixed -values;  E^  and  (dX/d0) 
are  constant  for  a  given  wavelength „  Therefore,  for  a  given 
wavelength  the  radiative  energy  emerging  from  the  exit  slit 
is  proportional  to  the  square  of  the  slit  width  setting.  Note 
that  Equation  (A-3)  was  derived  by  neglecting  diffraction. 
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APPENDIX  B 


DETERMINATION  OF  RESOLVED  BAND  WIDTH  (OR 
SPECTRAL  SLIT  WIDTH)  OF  THE  MONOCHROMATOR 

The  resolved  band  width,  AX,  of  a  monochromator  is  de¬ 
fined  as  half  of  the  total  wavelength  interval  falling  on 
the  exit  slit;  it  is  a  measure  of  the  wavelength  purity  of 
the  beam  emerging  from  the  exit  slit.  It  is  necessary  to 
know  the  resolved  band  width  as  a  function  of  slit  width  and 
wavelength.  This  function  will  aid  the  choosing  of  proper 
slit  widths  at  different  wavelengths  for  reflectivity  measure 
ments.  These  slit  width  settings  will  give  minimum  error  due 
to  the  effect  of  slit  opening  and  still  allow  enough  energy 
for  detection  (23). 

The  resolved  band  width  can  be  written  as 

AX  =  (AX)d+  (AX)p  (A-4) 

The  first  terra  (A  X)d  is  due  to  the  dispersion  of  the  prism 

and  the  second  term  (AX)  is  the  Rayleigh  term  and  results 

P 

from  the  diffraction  pattern  of  an  infinitely  narrow  slit. 
Both  terms  will  be  derived  in  detail  - 

Determination  of  (AX)^  -  The  spectral  range  of  radia¬ 
tion  going  through  the  prism  once  and  falling  on  the  exit 
slit  of  width  s  in  the  focal  plane  is  (11) 
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(A\)d  -  ■  •  A  0  -(|4)df 


which  can  be  written  as 

1 


(AX)„  = 


s> 

f 


(A-5) 


(A-6) 


In  the  case  of  a  prism  used  at  minimum  deviation  , 

in\2  / 


d0  _  2sin\2/ 


(A-7) 


dn  (l-n2sin2  (a/2))^ 
where  n  is  the  index  of  refraction  of  the  prism  at  wavelength 
A  and  oc  is  the  prism  apex  angle.  Therefore 

[l~n2sin2(a/2)jJs  ± 


(AA)d  2  sin  (a/2 f 


(dn7d XT  ’  f 


(A-8) 


For  a  single  pass  monochromator  the  radiation  passes 
through  the  prism  two  times  and  both  the  entrance  and  exit 

Vc  “fc 

slits  have  the  same  slit  width  ,  the  resolved  band  width 
due  to  the  dispersion  of  the  prism  should  be  half  of  the 
above ,  i .  e . 


i  h 


^l-n2sin2 (a/2) j 
(A  A)d=  4  sin  (a/2)  (dn/dA)  “  I 


(A-9) 


When  prisms  are  used  in  spectrometers,  they  are  always  set 
as  nearly  as  possible  at  minimum  deviation  because  otherwise 
any  slight  divergence  or  convergence  of  the  incident  light 
would  cause  astigmatism  in  the  image  (33) 

For:  (a)  The  monochromator  can  be  used  interchangeably  be~ 
tween  entrance  slit  and  exit  slit, 

(b)  For  a  given  sum  of  entrance  and  exit  slits,  the 
errors  or  corrections  due  to  wide  slits  are  least 
when  the  two  slits  are  equal  (23). 
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Determination  of  (AX)  -  The  equation  of  resolved  band 

ir 

width  due  to  diffraction  limitation  of  a  prism  is  well  known 
and  is  given  as  (33) 

(AX)P  =  b (dn/dX )  (A-10) 

where  b  is  the  effective  prism  base.  For  a  single  pass 
monochromator  the  radiation  passes  through  the  prism  two 
times,  which  will  give  two  times  better  resolution  as  men¬ 
tioned  above,  i.e. 

(A  X)p  =  2b(dn7dXT  (A-H) 

Actually  (AX)  will  depend  on  the  slit  width  as  well  (40) 

r' 

(AX)p  =  2b  (dn/dX )  (A-12) 

with  F(s)  a  function  of  slit  width.  Usually,  (AX)d  is 
much  larger  than  the  term  (A  X)^,  thus  F(s)  is  only  a  second 
order  correction.  For  simplicity  without  introducing  large 
error,  set  F(s)=l,  and  the  equation  of  (AX)  becomes 

r  o  2 

[l-n^sin  (a/2)J  g  x 

A  X  4  sin  (S’/ 2 )  ( dn/dXj  *  f  +  2b  (dn/dX)  (A-13) 

for  a  single  pass  monochromator.  Equation  ( A— 1 3 )  can  be 

written  as 

A  X  =  A  Xq^*  q—j-  +  A  Xp  (A-14) 

where  A  Xq  ^  is  the  resolved  band  width  due  to  dispersion 
for  0.1  mm  slit  width  and  s  is  the  slit  width  used  in  the 
unit  of  mm. 

In  Figure  8a  of  reference  (40),  A  Xg.i  and  A  Xp  as 
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functions  of  wavelength  are  given  for  the  double  pass  mono¬ 
chromator.  These  values  should  be  doubled  in  order  to  cal¬ 
culate  the  resolved  band  width  for  a  single  pass  monochroma¬ 
tor. 

By  using  values  of  A  Xq  ^  and  A  X^  taken  from  curves 
in  Figure  8a  of  reference  (40)  and  using  Equation  (A-14) , 
the  resolved  band  widths  of  the  P-E  Model  98  single  pass 
monochromator  are  calculated  for  wavelengths  from  two  to 
13  microns  and  with  slit  openings  from  0.1  mm  to  1.6  mm. 

The  results  are  listed  in  Table  A-l. 

Table  A-l.  Resolved  Band  Width  (in  microns) 
of  Single  Pass  Monochromator 


Wavelength  Slit  Opening  (mm) 


M 

0.1 

o 

• 

to 

0.4 

CO 

• 

o 

1.60 

2 

0.0450 

0.0850 

0.1650 

0.3250 

0.6450 

3 

0.0604 

0.1124 

0.2164 

0.4202 

0.8404 

4 

0.0672 

0.1042 

0.1982 

0.3862 

0.7622 

5 

0.0478 

0.0848 

0.1588 

0.3068 

0.6028 

6 

0.0420 

0.0730 

0.1350 

0.2590 

0.5070 

7 

0.0370 

0.0630 

0.1150 

0.2190 

0.4270 

8 

0.0338 

0.0568 

0.1028 

0.1948 

0.3788 

9 

0.0304 

0.0504 

0.0904 

0.1704 

0.3304 

10 

0.0288 

0.0478 

0.0858 

0.1618 

0.3138 

11 

0.G274 

0.0454 

0.0804 

0.1534 

0.2974 

12 

0.0264 

0.0438 

0.0786 

0.1482 

0.2874 

13 

0.0260 

0.0432 

0.0786 

0.1464 

0.2840 

APPENDIX  C 


THEORY  CONCERNING  FLUX  AVERAGING  SPHERE 


A  flux  averaging  sphere  is  a  sphere  with  diffusive  inner 
coating  and  used  over  the  detector  (see  Figure  A-2) .  The 
function  of  the  flux  averaging  sphere  is  to  spread  the  inci¬ 
dent  flux  uniformly  over  the  sensitive  area  of  the  detector, 
regardless  of  image  size,  shape,  or  intensity  distribution. 
One  can  thus  eliminate  the  problems  of  spatial  sensitivity, 
angular  sensitivity,  and  collection  of  all  of  the  incident 
flux  due  to  the  image  over  filling  the  detector  which  might 
be  caused  by  aberrations  and  any  slight  mis-alignment  of  the 
optical  components. 

Jacquez  and  Kuppenheim  (32)  developed  the  general  theory 
for  the  performance  of  an  integrating  sphere  expressed  in  the 
relation 


Bu  =  P 


b/S 

rj  dY  1 

s  / 


c/S 


=  (Prs)7? 


(A-15) 


where 


Bd  is  the  total  flux  passing  out  through  the  detector 
port  with  tiie  sample  aL  the  sample  aperture. 

P  is  the  total  radiant  incident  flux  into  the  entrance 


port 


o 


Incident  Flux 


Detector  Housing 

Figure  a-2.  Sketch  of  Flux  Averaging  Sphere. 


O 


7]  is  the  efficiency  of  the  sphere. 

r  is  the  reflectance  of  the  sample, 
s 

a  is  the  spherical  area  of  entrance  port. 

b  is  the  spherical  area  of  detector  port. 

c  is  the  spherical  area  of  sample  port. 

2 

S=4ttr  ,  R  is  the  radius  of  the  sphere. 

r^  is  the  reflectance  of  sphere  inner  coating. 

d=S-a-b-c 


An  averaging  sphere  has  only  one  entrance  port,  one 
detector  port,  and  no  sample  port,  hence  the  general  relation 
reduces  to  a  simpler  form  in  the  case  of  an  "averaging  sphere 


bd  =  p 


_ _ _ b/S  _ 

(1/r  ,)  -1  +  (a/s)  +  (b/s) 


(A-16) 


For  designing  a  sphere  to  serve  the  purpose  of  flux 
averaging,  the  sphere  should  fulfill  the  following  conditions 
1.  The  wall  of  the  sphere  should  be  coated  with  a  material 
which  will  give  high  reflectance  in  the  wavelength  range  con- 
cerned. 

2„  The  sphere  wall  should  be  a  diffusing  surface  if  a  con¬ 
stant  intensity  across  the  detector  port  is  to  be  expected. 

3.  The  sphere  size  should  be  small. 

4.  The  detector  port  should  be  as  large  as  practical  . 

5.  The  entrance  port  should  be  as  small  as  possible. 

Averaging  spheres  of  two  and  three  inch  diameters,  with 


Roto-blasting  inner  surface  and  sulfur  coating  (18)  were 
fabricated.  Both  entrance  and  detector  ports  were  of  the 
size  5/16  inch  x  1/2  inch. 
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With  the  above  given  dimensions,  the  total  flux  reach¬ 
ing  the  detector,  B^,  can  be  estimated  theoretically.  For 
simplicity  in  calculation  of  B^,  the  areas  of  entrance  and 
detector  ports  are  neglected  comparing  to  the  total  surface 
area  of  the  sphere.  Thus,  the  equation  of  B^  become 

bd  =  p  717^  tA'17) 

The  reflectance  of  the  coating  (sulfur),  r  ,  is  from  0.95 
to  0.75  for  the  wavelength  range  from  1.5M  to  7.0 M  (18). 

The  total  flux  reaching  the  detector,  B^,  is  calculated 
for  averaging  spheres  with  two  and  three  inch  diameter  at 
wavelengths  1.5/J  and  7.0 M.  The  results  of  B^  are  given  in 
percentages  of  the  incident  flux  to  the  entrance  port  (see 
Table  A-2). 


T<_ble  A-2.  Total  Flux  Reaching  The  Detector. 


M/0 

Total  Flux  Reaching 

the  Detector 

3  inch  sphere 

2  inch  sphere 

1.5 

11.0%  P 

25.0%,  P 

7.0 

1.7%  P 

3.7%  P 

Different  coating  thicknesses  and  detector  port  con¬ 
figurations  of  the  averaging  spheres  were  tried.  The  maxi¬ 
mum  radiant  flux  received  on  the  detector  was  about  one- 
fourth  of  the  value  oredicted  in  Table  A-2.  That  is,  at 
about  2M  the  radiant  flux  received  on  the  detector  with  the 
averaging  sphere  is  only  about  6%  of  the  radiant  flux 
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received  on  the  detector  without  the  averaging  sphere  and 
at  7.0/i  it  is  only  C<,9%. 


US 


VITA 


Jack  Jin,»  Cue  Hsia  is  a  citizen  of  the  Republic  of  China, 


After- high  school,  Mr.  Hsia  attended  National  Taiwan 
University,  and  received  the  B.S.M.E.  degree  in  June  1959. 
Upon  receipt  of  his  degree  from  National  Taiwan  University, 
he  served  as  a  ground  engineer  with  the  rank  of  second  lieu¬ 
tenant  in  the  Chinese  Air  Force  until  February  1961. 

in  September  1961,  Mr.  Hsia  entered  Purdue  University, 
.where  he  served  as  a  graduate  research  assistant  at  the 
■Thermophysical  Properties  Research  Center  and  was  awarded 
the  M.S.M.E.  degree  in  June  1964.  He  remained  at  Purdue 
University  to  study  towards  a  doctorate. 

Mr.  Hsia  is  a  student  member  of  the  optical  Society  of 
America  and  an  associate  member  of  the  American  Society  of 
Mechanical  Engineers. 

Mr.  Hsia  and  his  wife. 


